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NOMENCLATURE 


Symbol 

Definition 

Units 

A 

Area 

ft* 

a 

Denotes  average  value 

SpeciHc  heat 

Btu/lbm  °F 

d 

Thickness 

in. 

D 

Diameter 

ft 

DE 

Equivalent  diameter 

ft 

E 

Young’s  modulus 

Ibf/in.* 

f 

Subscript  referring  to  final  value 

Or 

Grashof  number 

h 

Heat  transfer  coefficient 

Btu/hr  ft*  ° 

H 

Height 

ft 

1 

Subscript  referring  to  ice 

k 

Thermal  conductivity 

Btu/hr  ft  °F 

L 

Latent  heat  of  fusion 

Btu/lbm 

m 

Subscript  referring  to  melting 

M 

Mass 

Ibm 

N, 

Heat  storage  parameter 

Nu 

Nusselt  number 

0 

Subscript  referring  to  initial  value 

Pr 

Prandtl  number 

Q 

Heat 

Btu 

Q 

Rate  of  heat  transfer 

Btu/hr 

R 

Radius 

ft 

Re 

Reynolds  number 

Subscript  referring  to  sink 

t 

Time 

hr 

T 

Temperature 

°F 

U 

Mean  flow  velocity 

ft/hr 

V 

Volume 

ft* 

w 

Subscript  referring  to  water 

W 

Mass  flow  rate 

Ibm/hr 

a 

Rate  of  change  of  temperature 

°F/hr 

e 

Strain 

in./in. 

0 

Stress 

Ibf/in.* 

P 

Density 

lbm/ft» 

0 

Dimensionless  time 

V 

Kinematic  viscosity 

ft* /hr 

CONVERSION  FACTORS:  U S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 

These  conversiDti  factors  include  all  the  significant  digits  given  in 
the  conversion  tables  in  the  ASTM  Metric  Practice  Guide  (F.  .^80), 
which  has  been  approved  by  the  Department  of  Defense.  Converted 
values  should  be  rounded  to  have  the  same  precision  as  the  original 


(sec  E 380). 

Multiply 

By 

To  obtain 

inch 

25.4* 

millimeter 

foot 

3.048* 

meter 

gallon/minute 

6.309020X  10'* 

meter’ /second 

pound-mass 

0.4535924 

kilogram 

pound-mass/foot^ 

16.01846 

kilogram/meter’ 

pound-mass/hour 

1.259938x10-^ 

kilogram/second 

pound-force /inch  ^ 

6.894757 

kilopascal 

Btu/hour 

0.2928751 

watt 

Btu/pound-mass  ’F 

4184.000 

joule/kilogram  kelvin 

Btu/hour  foot^  °F 

5.674466 

watt/meter’  kelvin 

degree  Fahrenheit 

f.^.Mr.^.-32)/1.8 

degree  Celsius 

EXPERIMENTAL  SCALING  STUDY  OF  AN  ANNULAR 
FLOW  ice-water  HEAT  SINK 

J.  Stubstad  and  W.  Quinn 

INTRODUCTION 

Tlie  investigation  of  efficient  methods  of  storing 
the  waste  heat  produced  by  the  generation  of  power 
in  hardened  underground  installations  has  resulted  in 
the  development  of  several  new  concepts  for  short- 
term heat  sinks.  Some  of  the  methods  explored  in- 
clude the  use  of  favorable  in-situ  conditions  to  pro- 
vide the  necessary  heat  storage  capacity.*  ’ ’’  ’’ 

For  example,  in  one  method,  steam  from  the  turbine 
is  condensed  on  the  surfaces  of  rock  tunnels  or  pas- 
sageways. Tlie  waste  heat  is  thus  “stored"  in  the  rock 
media  surrounding  the  tunnel.’  ’’ 

Another  concept  involves  the  use  of  an  ice-water 
heat  sink  which  is  relatively  independent  of  in-situ 
ground  conditions  (see  cover).  Waste  heat  is  stored 
by  the  conversion  of  ice  to  water  and  subsequent  heat- 
ing of  the  water.  Relative  to  other  types  of  liquid  reser- 
voir heat  sinks,  the  ice-water  sink  provides  a longer 
period  of  low  and  essentially  constant  temperature 
coolant  water  for  the  condenser.  This  results  in  high 
overall  thermal  efficiencies  for  the  power  plant,  par- 
ticularly during  the  more  critical  early  stages  of  sink 
usage.  In  addition,  the  large  amount  of  heat  required 
to  convert  ice  to  water  results  in  high  volumetric  effi- 
ciency. This  increased  storage  of  heat  per  unit  volume 
of  sink  reduces  construction  costs. 

Using  a 4-ft-diam  and  6-ft-high  scale  model,  CRREL 
initially  examined  the  feasibility  and  practicality  of 
the  ice-water  concept.*  A computer  program  was  de- 
veloped during  this  initial  study  to  numerically  simu- 
late the  temperature-time  history  of  the  model  sink. 

This  study  concluded  that  an  annular  flow  ice-water 
heat  sink  was  feasible  and  that  the  computer  program 


could  predict  the  behavior  of  the  scale  model  with  a 
reasonable  degree  of  accuracy.  A subsequent  study 
was  conducted,  using  the  same  model  but  examining 
different  ice  fomis,  specifically  ice  cylinders,  ice 
blocks  and  ice  cubes.”  It  concluded  that  the  use  of 
a solid  ice  cylinder  represented  the  best  overall  bal- 
ance between  a desired  low,  essentially  constant  out- 
let water  temperature  and  optimum  heat  storage 
capacity. 

A m jor  concern  recognized  during  these  initial 
studies  vas  the  potential  problem  in  scaling  the 
results  of  tMs  rather  small  model  to  a large  prototype. 
With  respect  n the  4-ft-diam  by  6-ft-higli  model,  a 
prototype  sink  (estimated  to  be  on  the  order  of  65  ft 
in  diameter  and  1 10  ft  high)  would  be  almost  5000 
times  larger  on  a volumetric  basis.  With  scaling  fac- 
tors of  this  magnitude  it  was  quite  possible  that  un- 
desirable characteristics  in  the  prototype  miglit  not 
be  observable  in  the  scale  model.  As  a result  another 
test  program  was  sponsored  by  the  Engineer  Power 
Group  (now  Facilities  Engineering  Support  Agency) 
to  investigate  possible  scaling  effects.  Space  limita- 
tions within  the  CRREL  facility  restricted  the  overall 
size  of  the  new  model  (Model  II)  to  6 ft  in  diameter 
by  10  ft  high,  yet  this  resulted  in  a volumetric  in- 
crease of  375%  over  the  previous  model.  Tlie  volu- 
metric ratio  between  a prototype  sink  and  the  new 
Model  II  apparatus  is  1300  to  1. 

In  addition  to  testing  scaling  effects,  the  Model  II 
test  series  also  investigated  some  of  the  practical  op- 
erational problems  that  had  been  encountered  during 
previous  test  programs,  for  example,  corrosion  pro- 
tection, economical  methods  of  creating  the  ice 
cylinder,  the  point  of  introduction  effect  of  inlet 
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Tins  report  represetits:  I ) the  experniiciital  results 
nsins;  the  laritei  Moilel  II  test  appaialus.  2)  a eonipari- 
son  ol  the  test  results  niili  those  predicted  hy  the  pre- 
vuuisly  developerl  conipnlei  ptoitiani.  and  2 ) a coni- 
pansrMi  ol  the  lesulls  with  tests  ciMidiicted  ni  earlier 
piopranis. 


DESCRII  nON  OK  EXPERIMENT 
Apparatus 

The  Model  II  expeiiiiienliil  apparatus  is  illustrated 
in  I igurcs  1 and  2.  It  represents  an  enlarged  version 
ol  the  model  used  in  the  previous  studies.  Due  to 
lesItictiiMis  in  laboratory  overhead  room,  the  height 
ol  the  new  model  (Model  II)  was  limited  to  10  It; 
a nominal  diamelet  ol  b H was  selected  to  provide 
a length  diamelet  ratio  similar  to  that  ol  a proto- 
type sink  IK)  It  high  and  bS  It  in  diameter.  Volu- 
nielrically  the  new  lest  apparatus  was  limes 
larget  than  the  pieviously  used  (i-H-higli,  4-lt-diani 
model  (Model  I). 

To  lree/,e  the  waiei  in  the  tank,  a tiichloroelhylene 
brine  at  a nominal  -7.2°T  temperature  was  circulated 
througli  SOO  It  ol  '4-iti.-diam  copper  tubing  wrapped 
around  the  sides  and  bottom  ol  the  lank.  Tlic  tank 
walls  were  insulated  with  .Tin.  higli-density  blanket- 
type  l iberglas  insulation  and  covered  with  a PVC 
vapor  barrier.  Two  layers  of  2-in.  sheets  of  polyure- 
thane were  used  to  insulate  the  bottom,  and  a paint- 
on  vapor  barrier  was  used  to  seal  the  bottom  and  cor- 
ners. Tlic  cover,  a composite  of  2x 4’s  and  lii-iti.  ply- 
wood. was  insulated  with  inlays  of  2-in.  polyurethane. 
Tlie  insulation  provided  an  excellent  barrier  to  envir- 
onmental iiinueiices,  yielding  an  overall  heat  transfer 
coclTicieni  ofonly  0.  IKio  0.20  Btu/hr  ft’  °F  be- 
tween the  oulsidc  and  inside  surfaces.  Tlie  heat  trans- 
fer coelTicicnt  for  the  insulated  apparatus  was  deter- 
mined by  filling  the  tank  with  hot  water  and  monitor- 
ing the  decrease  in  hulk  water  temperature  during  a 
24-hr  period.  Combined  with  air  temperature  read- 
ings during  this  penod,  the  decrease  in  sensible  heat 
was  converted  into  a mean  heat  transfer  coefficient. 


To  prevent  possible  development  ol  the  htgli 
stresses  associated  with  the  enlrapmeiit  ol  hi|uid 
walei  hy  foimnig  ice.  a system  was  installed  to  bub- 
ble air  through  the  center  ot  the  tank  and  acioss  the 
surfaces  ol  the  viewports.  The  air  sticam  provided  a 
coiiiinuons  free  path  to  the  surface  and  iheiehy  pre- 
vented the  entrapment  of  liquid  water  by  the  ice. 

During  the  previous  model  iesis  dilTicullies  were 
experienced  in  providing  corrosion  proleclioii,  due  to 
the  extreme  variations  in  tank  wall  temperature  ol 
-70°r  to  ■l■120°E.  Large  stresses  produced  by  the  dil- 
feieiit  coeliicients  of  thermal  expansion  I'oi  the  sur- 
face coaling  and  the  steel  wall  resulted  m the  scpaia- 
tion  of  the  linisli  from  the  wall  and  the  loss  ol  corro- 
SKMi  protection.  Discussions  with  manulacturers  ol 
various  coirosion-resistant  paints  indicated  that  no 
staiidaid  paint  was  maniifacturcd  which  could  with- 
stand both  the  extreme  low  and  high  lemperaiuies 
and  high  humidity  conditions  of  the  model  sink.  Tlie 
sohitiiMi  was  fouiul  througli  the  application  ol  a cold 
galvani/ing  compound  to  the  interior  surfaces.  Tliis 
product,  manulacturcd  by  the  ZRC  Corporation,  is  a 
mixture  of  approximately  /.iiic  and  .‘5%  binder, 
which  provirles  a sirrface  rinish  having  the  same  basic 
characteristics  as  hot  dip  galvani/ing.  fii  the  six  ex- 
[xtrimental  tests  conducted  during  the  Model  II  pro- 
gram. the  (Mily  failures  of  this  coaling  (Kcurred  m 
areas  where  the  steel  surface  had  not  been  sulTicienily 
cleaned  prior  to  coating.  A coating  of  this  type  would 
be  recommended  for  a prototype  steel  tank. 

As  with  the  Model  I apparatus,  two  viewports  were 
installed  in  the  tank  wall.  Observations  of  the  melt- 
ing geometry  were  made  by  using  the  viewports  and 
an  access  opeiimg  in  the  tank  cover.  A probe  assem- 
bly, inserted  througli  the  cover,  measured  the  chang- 
ing shape  of  the  ice  cylinder  during  a test,  Tlie  probe 
consisted  of  a stainless  steel  rod  and  a mounting  hlrK  k 
to  l(Kate  it  on  the  cover;  the  rod  was  lowered  into 
the  tank  until  it  touched  the  ice  cylindei.  By  meas- 
uring the  depiii  of  the  ice  from  the  tank  cover,  the 
heiglit  of  the  ice  at  that  point  was  determined.  Tlie 
probe  was  then  moved  outward  along  one  of  the 
three  radial  slots  in  the  cover.  One  slot  was  located 
directly  above  one  of  the  water  inlet  pipes  to  pemiit 
examination  of  the  scouring  pattern.  The  two  other 
slots  were  located  on  a common  diametei  but  be- 
tween two  outlet  pipes  to  avoid  as  much  of  the 
scouring  influence  as  possible. 


I 'iintrc  .1  Water  level  inJieating  glass  and  viewport  ( eoveredl 


i'igure  4.  Inlet  header. 


Figure  5.  Thermoeouple  arrangement 


Figure  6.  Interior  view  of  the  heat  sink  apparatus. 


Figure  7.  lee  hloek  loading  of  the  Model  II  sink. 
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('iKiiiges  in  water  level  were  monitored  using  a 
water  level  indieating  luhe  (f'lg.  .1).  Tlicse  changes 
were  then  related  to  the  volumetric  contraction  as- 
siKiated  with  conversion  of  the  relatively  low  density 
ice  to  water,  as  an  indication  of  the  gross  melting  rate. 

To  permit  comparison  with  the  Model  I test  series, 
acential  coolant  water  outlet  and  a six-point  inlet 
header  were  again  used  in  this  series  (Fig.  4).  Althougli 
mote  efficient  methods  of  introduction  and  removal 
of  the  coolant  water  liad  been  proposed,  it  was  desired 
to  replicate  the  pievious  test  conditions  as  much  as 
possible. 

Temperature  vaiiations  within  the  sink  were  meas- 
ured by  using  three  strings  of  thermocouples,  one 
placed  along  the  vertical  centerline,  and  the  others 
along  vertical  lines  at  the  ‘/i  and  X radial  position 
(F-ig.  5).  Tlie  thermocouples  were  fixed  in  position  by 
attaching  them  to  ropes  tied  to  anchor  points  at  the 
top  and  bottom  of  the  tank.  Tlic  thermocouples  were 
monitored  with  a Leeds  and  Northrup  chart  recorder, 
with  an  extra  thermocouple  in  an  ice  bath  for  calibra- 
tion purposes.  Inlet  and  outlet  water  temperatures 
were  measured  using  thermocouples,  two  in  the  inlet 
header  and  two  in  the  outlet  pipe;  a Kaye  data  logger 
was  used  to  take  readings  at  20-min  intervals. 

A 24-kW  immersion  healer  was  used  to  simulate 
the  power  plant  condenser;  control  of  the  rale  of  heat- 
ing was  provided  by  a Variac  Iransfonner.  A Vi-hp 
variable-speed  FX'  motor  with  a rotary  displacement 
pump  was  used  to  circulate  water  througli  the  heater 
and  back  into  the  sink.  Flow  rates  were  monitored 
with  either  a 0.5  to  6.0  gpni  or  1 .72  to  1 7.2  gptii 
rotameter,  the  latter  being  used  only  when  test  How 
rates  exceeded  5.0  gpni.  Tlie  heater  and  pump  mod- 
ule is  illustrated  in  Figure  2b. 

Oice  the  sink  had  been  completely  fro/.en.  an  ini- 
tial annulus  was  created  by  circulating  hot  brine 
througli  the  refrigeration  tubing.  A 6-kW  immersion 
healer  was  used  to  heal  the  brine. 

Tlie  use  of  hot  bnne.  rather  than  electric  heating 
tapes  which  were  employed  with  tlie  previous  appara- 
tus. was  chosen  because  this  type  of  system  provides 
maximum  utih/ation  of  the  refrigeration  lines  and  is 
believed  to  better  simulate  probable  heating  and  cool- 
ing schemes  for  a prototype  sink. 

Dunng  a test,  motion  of  the  ice  cylinder  was  pre- 
vented by  an  anchor  system  consisting  of  several 
boards  laslened  to  a rope  positioned  along  the  verti- 
cal centerline  of  the  lank.  WFiilc  some  rotational  mo- 


tion from  uneven  melting  did  occur,  the  gross  overall 
movement  was  insignificant.  Tlie  ice  anchor  system 
and  iherniocouples  are  shown  in  Figiiie  6. 

Test  procedure 

Two  techniques  were  employed  to  form  the  ice 
cylinder.  The  primary  method  involved  tilling  the 
test  lank  with  water  and  freezing  inward  from  the 
lank  wall.  Approximately  six  days  were  required  to 
completely  freeze  all  the  water.  Two  tests  were  con- 
ducted with  a second  method,  the  loading  of  a mixture 
of  ice  and  water  prior  to  the  start  of  freezing.  Tliese 
tests,  employing  an  initial  ice  load  of  6000  to  8(XKI  lb 
in  50-lb  blocks  and  a water  load  of  6(X)0  to  SCKX)  lb. 
reduced  the  freezing  period  to  about  five  days.  A 
view  of  the  heat  sink  during  the  block  loading  opera- 
tion is  shown  in  Figure  7.  Tlie  cost  effectiveness  of 
this  priKedure  is  somewhat  obscured  by  the  availa- 
bility of  a large  capacity  on-line  refrigeration  system 
at  CRRFvL.  but  even  an  initial  starting  mass  of  50'^ 
ice  does  not  significantly  affect  the  freezing  period. 

All  explanation  is  that  the  longest  part  of  the  freez- 
ing period  is  associated  with  freezing  the  center  of  the 
ice  cylinder,  and  an  initial  load  of  ice  does  not  increase 
the  rate  of  heat  conduction  througli  the  ice.  Of  course 
this  procedure  does  result  in  a substantial  decrease  in 
the  on-site  refrigeration  load.  (The  possibility  of  gen- 
erating higlier  packing  ratios,  i.e.  initial  ice  to  total 
mass  ratios,  was  investigated  in  a previous  study,'* 
and  it  was  concluded  that  for  common  ice  shapes, 
such  as  cubes  or  blocks,  a maximum  of  about  61% 
ice  was  the  best  that  could  be  obtained.)  A detailed 
analysis  of  the  freezing  process  is  given  in  Appendix 
D. 

Once  freezing  had  been  completed,  the  sink  was 
allowed  to  stand  at  room  temperature  for  a period  of 
one  to  two  days.  A temperature  gradient  of-70°F 
at  the  tank  wall  to  32°F  at  the  center  of  the  tank  nor- 
mally existed  at  the  moment  when  complete  freezing 
had  occurred;  this  standby  period  allowed  tempera- 
ture equilibration  to  exteur.  Wlien  the  tank  wall  had 
reached  a temperature  of  approximately  -10°F,  the 
brine  heater  was  activated  and  the  melting  of  an  ini- 
tial annulus  commenced.  Since  the  hot  brine  llowed 
around  the  tank  througli  one  continuous  loop,  it  was 
necessary  to  reverse  the  direction  of  brine  flow  sev- 
eral times  in  order  to  melt  a uniform  annulus.  Tliis 
technique  was  only  partially  successful;  the  initial 
annulus  nomially  varied  from  I to  2 in.  at  the  top 
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to  '-i  to  '4  111.  at  the  hollom  corner.  In  Inlnrc  experi- 
ments the  brine  system  shoiiM  be  composeil  of  several 
circuits  connecleil  m a parallel  manner  to  ensuie  de- 
velopment ol  a more  unilorm  annulus. 

Alter  the  initial  annulus  had  been  lormed  aiul  the 
average  ice  temperature  was  between  and  .^2°l  . 
water  was  added  to  completely  submerge  the  ice  evhn- 
dei.  Bv  covering  the  ice  bh>ck  with  water  during  the 
lest,  changes  in  the  water  level  could  be  relaieil  to  the 
amount  of  ice  melted  at  anv  lime.  Alter  a short  peri- 
od of  lime  to  allow  this  added  water  to  reach  an  et|ui- 
libnuin  temperature,  the  lest  was  begun. 

IXiring  the  test  the  coolant  water  How  rale  and  the 
temperature  difference  between  the  sink  inlet  and 
outlet  water  were  monitored  to  check  the  rale  of  heal 
rejection.  Readings  of  the  change  in  water  level  and 
shape  of  the  ice  were  taken  at  random  intervals.  Once 
the  ice  was  completely  melted  and  the  outlet  water 
temperature  was  at  least  bO°F.  the  lest  was  slopped 
and  a measurement  was  made  rrf  the  bulk  water  temp- 
erature. From  this  and  the  initial  ice  mass  and  temp- 
erature. a heat  balance  was  computed  as  a check  on 
the  nominal  heal  rejection  rate. 

EXPERIMENTAL  RESULTS 

Suminary  of  experimental  tests 

Six  expenmental  tests  were  conducted  with  a nom- 
inal heat  rejection  rate  of  42,01.^  Btu/hr.  Coolant 
water  flow  rates  varied  from  2.50  to  13.50  gpm.  Both 
the  heat  rejection  rate  and  coolant  water  flow  rates 
were  selected  by  using  the  scaling  relationships  de- 
veloped by  Brown  and  Quinn.’  Table  I summarixes 
the  tests. 

By  measurement  of  the  ice  cylinder  at  the  start  of 
a test  and  the  total  change  in  sink  volume  at  the  com- 
pletion of  melting,  the  initial  ice  mass  was  determined. 
In  the  six  tests  the  initial  ice  mass  varied  from  a low 
of  12,250  Ibm  in  test  2 to  a higli  of  13,420  Ibm  in 
test  5,  the  ice  mass  averaging  12,970  Ibm.  The  total 
mass  (both  liquid  and  solid  states)  contained  in  each 
test  varied  from  14,150  Ibm  in  test  4 to  14,370  Ibm 
in  test  3,  with  an  average  of  14,280  Ibm.  Of  critical 
importance  to  the  volumetric  efficiency  of  an  ice- 
water  heat  sink  is  the  percentage  of  total  mass  in  the 
form  of  ice;  in  this  test  series  it  was  found  that  the 
ratio  varied  from  a low  of  85.7%  to  a higli  of  94.0% 
with  an  average  load  factor  of  90.8%. 


Tliroe  separate  procedures  were  emploved  lor  com- 
(larison  loeslabltsh  the  average  treat  reieclron  rale  dur- 
ing the  tests 

/‘roccdiirc  /.  With  the  heal  sink  inlet  and  outlet 
coolant  water  lemjx’raliires  and  the  Ilow  rate  known 
at  any  particular  instant,  the  heat  reiection  late  (J  can 
be  determined  from 

= ■ ^-ui*  (Btu/hr) 

where  h'  = mass  How  late  (Ibrn/hr) 

= spectfic  heat  ol  water  (Btu/lbrn  °F  ) 

^ui  ■ temperature  (°F) 

= outlet  temperature  (°F). 

To  determine  the  average  value  of  the  heal  rejection 
rale  Q.^  during  the  test,  the  following  lormiila  can  be 
used  lor  “n"  time  intervals: 

= E <^..-^our)  <') 

where  the  ma.ss  How  rale  Iv'and  sjKCilic  heal  r p have 
been  assumed  to  be  constant.  In  practice,  readings  were 
taken  at  20-min  intervals  so  that  the  average  value  is 
based  on  150  to  2(K)  separate  temperature  readings. 

Proccdurv  2.  At  the  completion  ol  a lest  a heat  bal- 
ance was  conducted  to  examine  the  reliability  of  the 
previous  method.  In  this  method  three  heal  storage 
components  must  be  considered:  sensible  heating  of 
the  ice  (>j,  latent  heat  , and  sensible  heating  of  the 
water  Q^.  Tlie  initial  ice  temperature  was  estimated  by 
averaging  readings  from  the  thermiK'ouples  in  the  ice 
(12  thermocouples  in  all)  and  assuming  that  this  aver- 
age value  was  the  mean  initial  ice  temperature  Tj.  Thus 
the  sensible  heat  storage  capacity  of  the  ice  is 

Q,=M,c„i(32-T0  (2a) 

where  is  the  initial  ice  mass  (Ibm),  and  Cpj  is  the  spe- 
cific heat  of  ice  (Btu/lbm  °F).  Tlie  latent  heat  storage 
capacity  for  the  melted  ice  is 

(2b) 

where  1.  is  the  latent  heat  of  fusion  (144  Btu/lbm). 
Finally  the  sensible  heat  storage  capacity  of  the  water 
is 
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Table  I.  Summary  of  tests;  heat  sink  Model  II. 

Tank  si^e;  (%  ft  diam.  10  ft  high; 

Nttminal  heat  rejection  rale:  42.013  Btu/hr. 


i\ominai 

Initial 

Tidal  mass 

Ratiif  ice! 

Test 

flow  rate 

ice  mass 

ice  water 

total  mass 

no. 

ixpmt 

(Ibrn) 

flhm} 

(%) 

1 A 

4.73 

13.210 

14,300 

92.4 

IB 

4.73 

13.290 

14.300 

92.9 

2 

2.50 

12,250 

14.300 

85.7 

3 

10.00 

12,750 

14,370 

88.7 

4 

7.50 

12.900 

14.150 

91.2 

5 

13.50 

13.420 

14.270 

94.0 

Average  values 

12,970 

14.280 

90.8 

Table  II.  Average  heat  rejection  rates.  Model  II. 

Nominal  heal  rejection  rate: 

: 42,013  Blu/hr. 

Nimnnal 

Q 

Q 

Q 

7>-v/ 

flow  rate 

Procedure  1 

Procedure  2 

Procedure  J 

no 

Ifpml 

(Btufhrj 

(Bti4/hr) 

(Btufhr) 

1 A 

4.73 

41,305 

38,570 

38.600 

IB 

4.73 

40,150 

35.050 

39,340 

2 

2.50 

41.140 

40,720 

29.060* 

3 

10.00 

41.410 

38,530 

36,240 

4 

7.50 

42.090 

35.920 

38,090 

5 

13.50 

40,260 

t 

34,790 

Average 

4 1 ,060 

37,760 

37.410* 

Deviation  from 

-2.3% 

-10.1% 

1 i.0% 

nominal  rate 


* Value  of  test  2 not  included  in  averatte  (see  text), 
t Not  available;  test  ended  prematurely  due  to  equipment  failure. 


whetc  Tf  = final  water  temperature  at  the  etui  of  the 
test  rn 

T„  = initial  water  temperature  (°F) 

- mass  of  water  in  the  annulus  at  the  start  of 
the  test  at  an  initial  temperature  of  Tq 
(Ibin). 

Tlie  average  heat  rejetlion  rale  by  the  heat  balance 
nKMhod  IS  determined  from 


(2d) 


where  l,  = time  to  reach  final  water  temperature  T, 

(hr). 

PniceJure  3.  This  rnclhcKl  considered  the  sensible 
heating  of  the  water  within  the  sink  after  all  ice  was 
completely  melted.  It  is  assumed  that  the  rate  of 


temperature  change  of  the  entire  sink  was  the  same  as 
the  rate  of  change  of  the  outlet  coolant  water;  that  is 

dt  dt 

where  is  the  bulk  “mixing  cup”  temperature  of  the 
sink. 

Such  an  assumption  represents  a good  approxima- 
tion only  when  there  is  a higli  degree  of  mixing  with- 
in the  sink  and  is  not  valid  when  density-induced 
stratificalion  occurs.  A least  squares  fit  is  made  to 
the  outlet  water  temperature  curve.  Tlie  slope  of  this 
line  isc/r,^,/r/r  so  that  the  siraiglit  line  heat  rejection 
rale  4?,  is 


pw  dt 


7 


whiMC  IS  iIk’  U)t;il  mass  of  watci  coiilainoil  within 
ilic  sink. 

Tlie  values  nf  the  average  heal  rejeclum  rale  using 
these  three  nielluHls  have  been  tahiilateJ  in  Table  II 
with  the  I'vcrall  average  lot  all  tests  ciiinbinetl.  As 
noleil.  the  third  method  eannol  be  used  when  den- 
siiy-relaied  siraiitleation  elTecis  predominate,  siieh  as 
in  lesi  2.  For  this  test  a heal  reieetion  rale  based  on 
the  slope  of  the  outlet  water  temperature  curve  is 
meaningless  and  has  not  been  included  in  the  overall 
average. 

A close  correlation  exists  between  the  overall  aver- 
age values  derived  by  Procedures  2 and  ,V  It  is  fell  that 
the  actual  experimental  heat  rejection  rales  lie  between 
the  maximum  and  minimum  limits  established  by  these 
three  methods.  Thus  the  experimental  heat  rejection 
rales  were  between  2 and  I I9f  less  than  the  nominal 
value  of  42.01. ^ Btu/hr. 

In  lire  three  procedures  used  to  calculate  the  average 
heat  refection  rale,  heal  transfer  to  the  surroundings 
was  assumed  to  be  negligible  and  was  not  included.  The 
validity  of  this  assumption  can  be  shown  by  examining 
the  results  of  lest  .V  During  this  test  the  average  air 
temperature  was  63.3‘’F.  A thermocouple  mounted  on 
the  inside  tank  wall  at  one-half  its  heiglit  indicated  an 
average  tank  wall  temperature  of  53.8°F  for  the  test. 
Assuming  that  this  temperature  is  indicative  of  the 
mean  lank  wall  temperature,  the  resulting  average  temp- 
erature difference  of  9.5‘’F  across  the  tank  insulation 
results  in  a heat  gain  of  only  500  Btu/hr.  Such  a rate 
is  approximately  1%  of  the  nominal  applied  heat  load 
and  therefore  is  insignificant  in  comparison  to  the  ap- 
plied heating  rate. 

Influence  of  coolant  water  flow  rates 

Over  the  range  of  flow  rates  tested  it  was  found  that 
a well-defined  transition  rxeurs  in  the  pattern  of  water 
flow  through  the  sink.  At  the  lowest  flow  rates  distinct 
water  layers  were  formed,  resulting  from  differences  in 
density  of  the  inlet  water  and  the  water  within  the  sink. 
Tlie  warmer  inlet  water  formed  a layer  on  top  of  the 
colder  water  within  the  sink,  with  a rather  distinct 
boundary  between  them.  Tltis  boundary  proceeded 
very  slowly  down  the  tank  and  produced  an  almost 
top-t(^bottom  melting  pattern.  Tliis  pattern,  illustrated 
in  Figure  Ka,  shows  the  ice  profile  fur  test  2 conducted 
at  a coolant  water  flow  rate  of  2.50  gpm.  It  should  be 
noted  that  the  horizontal  upper  ice  surface  indicates  an 
almost  constant  rate  of  melting  across  the  entire  upper 


suil.icc.  even  though  the  healcil  iiilcl  walei  was  intro- 
duced at  the  iwo-thiids  radial  positiou 

This  phenoiiieiioii  is  also  apparent  alter  all  lire  ice 
has  melleil.  as  can  be  seen  m the  lime-lcmperalure  his- 
loiy  ol  the  outlet  water  lemperaluie  ( Appendix  l iguie 
A2 1.  \ole  ih.il  when  incllmg  is  completed  the  oiillel 
w.ilei  leinpeialuie  rapidly  iiicieases  amt  then  slabih/es 
al  a coiisiaiil  lemperaliire  which  erpials  that  id  die  in- 
let  walei  piioi  lo  the  sudden  inciease.  Il  then  lemaiiis 
al  ihis  lempeialiire  loi  appioximalely  10  hours  when 
It  again  shows  a sharp  mciease  B\  assuming  a umlomi 
velocity  prolile.  the  flow  rale  ol  2.50  gpm  ihroiigh  a 
b fl-diam  cyhndei  would  result  m a net  dowiiw.iid 
velocity  ol  0.142  in. /min.  L'ndei  conditions  wheie 
there  is  absolutely  no  mixing  beiween  l.iyeis.  appioxi 
maiely  I 1.5  hours  would  be  iec|uiied  loi  this  theo- 
retical layci  to  traverse  the  9b  in.  Iioni  the  walei  sur- 
face to  the  outlet.  The  actual  peiiod  ol  tins  tempera- 
ture plateau  was  10.5  hours  VViih  a modest  allowance 
for  some  mixing  of  the  water  along  the  boundary  . 
which  would  reduce  the  lime  rei|uiied  lor  the  holtei 
inlet  water  to  influence  the  outlet  water  lemperaluie. 
il  is  found  that  Ihc  expenmeulal  lesults  agree  lathei 
well  with  such  a flow  piixess. 

Tims,  al  low  flow  lales  dining  this  posl-meltmg 
stable  icmpeialiiie  period,  the  outlet  water  tempeia- 
ture  IS  essenlialK  nnallecled  by  the  addition  of 
heated  ml  ■!  water,  dial  is.  die  outlet  water  tempera- 
ture will  be  the  same  wheihei  the  coiulenser  water  is 
leliiineil  to  the  sink  oi  waster!  Tlie  duialion  of  this 
plateau  peiiod  can  he  expressed  simply  as 

M 

/ = 0‘)_'  (hr)  (4| 

li' 

where  B'ls  the  nia.ss  flow  rate  id  water  m Ibni/hi.  The 
factor  of  0.9  is  used  lo  account  loi  mixing elfects  al 
the  boundary. 

It  should  be  emphasized  that  this  formula  applies 
only  lo  low  coolant  water  flow  rates.  2.5  gpm  and 
less  in  the  Model  II  sink,  and  only  under  conditions 
of  constant  coolant  flow  and  heat  rejection  rales. 

If  we  attempt  lo  apply  the  same  approach  to  the 
tests  having  a flow  rate  of  4.73  gpm.  the  foiiiiula  pie- 
dicts  that  fur  the  nominal  heat  rejection  rate  of 
42,013  Btu/hr  there  should  be  a post-melting  stable 
temperature  period  lasting  S.S  hours.  A temperature 
plateau  of  this  type  did  not  occur;  instead,  changes 
of  slope  in  the  outlet  water  temperature-time  curve 
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iIi'M'Ii'IhhI  ,il  live-  lo  si\  lixtii  mlciAals.  indicaling 
lli.il  lliiMi-  w.iv  sDiiK'  lcinlcni.y  lor  sltalitUation  (I'lg. 

Al  ) TIiom'  poiunK  also  mdicalc  a miMiigol  warmer 
miel  walei  willi  ihe  crrirlet  sink  walei.  It  would  ap- 
|vai  that  ilensiis  iinliKcd  slialifKalion  is  not  the 
ilommaiit  Ilow  pioeess  tlirougli  the  sink  at  this  flow 
late,  hut  ijihei  iheie  is  prohably  an  almost  equal  bal- 
aiKe  between  oseiall  mixing  and  sliatiflealion.  Tims 
the  evidenee  ol  thermal  stratilieation  prevalent  in  the 
’ 50  epm  lest  is  shaiply  rerluced  al  higliei  flow  rates, 
and  the  iiansition  Irom  slralified  lo  mixed  flow  ck- 
kUis  m the  neighhorhood  of  4 75  gpm  tor  Ihe  Model 
II  apparatus 

l■ollowmg  the  same  reasoning  it  would  be  ex|ieeied 
that  the  UK  tease  m flow  rate  rrom  4 7.5  to  7.50  gpm 
wtnild  result  m Ihe  total  disappearance  ol'density- 
lelated  effects,  with  mixed  flow  being  the  doiitinani 
mode.  As  shown  in  Appendix  l igure  A.5  Ihe  outlet 
walei  temperature  curve  lor  Ihe  7.50-gpm  lest  does  be- 
come a slraiglil  line  after  completion  of  melting.  Ihe 
mcrcasing  linearity  resulting  from  Ihe  mixing  of  warm- 
er inlet  water  with  the  cohlei  water  within  the  sink. 
(The  terms  nanncr  uikI  calder  are  relative,  since  in  this 
lest  there  is  only  an  1 1°F  difference  in  temperature  be- 
tween Ihe  outlet  and  inlet  water.) 

This  behavior  is  also  apparent  at  Ihe  two  higliest 
flow  rales  used:  10.00  and  1.5.50  gpm.  For  flow  rates 
greater  than  7.50  gpm  in  the  Model  II  apparatus,  the 
sink  exhibits  a mixed  flow  condition  after  completion 
ol  melting.  As  in  the  low  flow  rate  tests,  this  leads  lo 
a lather  simple  lormula  for  predicting  the  post-melting 
lempeialure-iime  behavior  the  heating  of  a well- 
mixed  large  mass  of  water.  In  parametric  form  this  is 
expressed  by 

J', ,,„(/)  = u(/-r,„)+T^  (°F)  (5a) 

where  r = running  lime  variable  (hr) 

= lime  al  which  melting  was  completed 
7"^  = outlet  water  temperature  at  Ihe  completion 
of  melting  (°F) 

a = the  rale  of  temperature  change  of  the  sink 
water  (°F/hr) 

where 

a = — (5b) 

M.c„ 


Tins  equation  may  be  used  for  either  a constant  or 
variable  laie  of  heal  rejection;  however  , in  Ihe  latter 
case  a will  also  become  a vanabic  and  the  graph  of  the 
outlet  temperature  will  be  a curved  rather  than  slraiglil 
line. 

As  Slated  previously,  further  increases  in  flow  rate 
(i.e.  lo  10.00  and  13,50  gpm)  did  not  result  in  any 
changes  in  the  shape  of  Ihe  outlet  water  temperature 
curve  but  only  affected  Ihe  mean  outlet  water  tempera- 
ture. 

In  coniunclion  with  the  change  from  stratified  lo 
mixed  flow  at  4.73  gpm  there  was  also  a significant 
change  in  the  melting  pattern.  In  the  2.50-gpm  flow 
rale  lest,  ice  melting  proceeded  vertically  from  Ihe  lop 
downward  with  a uniform  rate  of  melting  across  the 
top.  Tlie  melting  pattern  for  the  lO.OO-gpm  test  was 
predominantly  radial.  A.s  might  be  anticipated  frrnn 
the  stratified  and  mixed  flow  conditions  which  de- 
veloped, the  melting  geometry  of  the  4.73-gpm  lest 
shows  both  vertical  and  radial  melting  patterns.  This 
would  indicate  that  stratified  flow  produces  a verti- 
cal melting  (top-downward)  pattern,  whereas  mixed 
flow  produces  radial  melting. 

Tltese  flow  rate  effects  can  also  be  illustrated  by 
examining  the  melting  pattern  of  the  ice  directly  be- 
neath one  of  the  inlet  pipes  (Fig.  8).  In  all  tests  the 
ice  was  completely  submerged.  In  the  2.50-gpm  test 
the  inlet  water  tends  to  create  a smooth  surface.  Tltis 
indicates  that  the  inlet  water  does  not  come  in  direct 
contact  with  the  ice,  except  during  the  first  few  hours, 
but  is  separated  from  it  by  the  water  above  the  ice  sur- 
face. This  is  not  true  for  the  4.73-gpm  flow  rate,  as 
shown  in  Figure  8b.  From  the  start  of  the  test  the 
inlet  water  produces  a high  rate  of  melting  near  the 
inlet,  as  evidenced  by  the  ice  scour  pattern.  Tltis  ef- 
fect has  a decreasing  influence;  by  the  twenty-second 
hour  the  melting  pattern  begins  to  resemble  the  geo- 
metry of  the  2.50-gpm  test.  Tlie  scouring  effect  be- 
comes increasingly  important  as  Ilow  rate  is  increased, 
so  that  for  the  lO.OO-gpm  test,  shown  in  Figure  8c,  it 
predominates.  After  just  eiglit  hours,  melting  is  al- 
most completely  radial  with  very  little  change  in  Ihe 
ice  height.  This  cutting  effect  of  the  inlet  water  is 
well  illustrated  in  Figure  9,  which  shows  the  ice  sur- 
face during  Ihe  lO.OO-gpm  test.  Note  Ihe  almost  ver- 
tical slice  in  the  ice  surface  in  the  vicinity  of  the  inlet. 
This  indicates  that  at  this  flow  rate  the  inlet  water 
causes  melting  along  a significant  length  of  the  ice 
cylinder. 
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Aiiollict  means  to  mvesiijtale  llic  Inllucncc  of  the 
eoolatit  water  Ilow  rate  ts  to  corttpate  outlet  teiiipeta- 
luie  time  histones  loi  the  tests,  as  shown  iti  (■'igure 
10  During  the  melting  stage  the  test  results  itiay  be 
partitioneil  Into  two  general  categories:  one  lot  the 
10  (K)-  atid  I .V50-gptn  tests  and  the  other  I'or  the 
2.50-,  4.7.1-  and  7.50-gptn  tests.  I'or  the  two  higitet 
flow  rate  tests,  there  is  a rather  smooth  rate  ol' 
chatige  in  outlet  water  temperature  starting  at  about 


4I°I' atid  hneatly  increasitig  at  aconstatit  rate  of  ap- 
proximately '/«  degrec/hr.  Kor  the  three  low  flow 
rate  tests,  there  are  fluctuations  in  the  outlet  water 
temperature  on  the  order  of  I to  .1  degrees  b'ahren- 
heit  within  an  overall  operating  range  of  36°  to  40°F. 
Since  the  maximum  density  of  water  occurs  at  3‘).2°F, 
it  is  possible  that  these  observed  fluctuations  are  the 
result  of  an  unstable  flow  condition  due  to  density  in- 
version effects.  This  would  indicate  that,  up  to  7..3-gpm 
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flow  rales,  the  outlet  water  is  lulluemeil  by  density  result  in  a lower  outlet  water  icnipeialure;  tiowevei. 

inversion,  the  efTects  ol  which  are  apparently  insig-  the  signiricanlly  lower  average  oiiilei  waiei  tempera- 

niticanl  at  the  higher  flow  rales.  lure  at  this  flow  rate  was  ipiiie  piohahly  inlliieiiceil 

Tlie  average  outlet  water  temperature  I'oi  ilie  entire  more  by  density-related  eirecis.  Tlie  possible  occiir- 

tnelting  period  vs  the  nominal  water  flow  rale  has  been  rence  of  a low  heal  rejection  rale  and  its  influence  are 

plotted  in  Figure  I la  and  tabulated  in  Table  III  It  can  discussed  in  later  sections. 

be  seen  that  except  for  the  7.50-gpm  lest  the  average  Average  outlet  water  temperatures  during  the  meli- 

oullcl  water  leinperalure  lies  along  a rather  smooth  "'g  period  vs  coolant  water  flow  rales  lor  the  Model  I 

curve  which  reaches  a maximum  value  of  4.3. S°F  foi  study  arc  shown  in  Figure  I lb,  Tliese  results  indicate 

flow  rales  of  about  10.00  gpni  and  higlier.  Tlie  dip  a maximum  average  outici  waier  temperature  con- 

which  (Kcurs  at  7.50  gpm  could  have  resulted  from  dition  developed  at  the  intermediate  flow  rates.  In  con- 

either  density-related  effects  or  from  a low  heat  rejec-  Itasi  the  results  ol  the  Model  II  study  indicate  a mini- 

tion  rale  during  the  test.  Wliile  the  possible  ixrcurrence  nitim  temperature  at  an  intermediate  flow  rale.  Tlie  ex- 

^ of  a low  heal  rejection  rale  is  not  indicated  in  Table  II.  isicnee  ol  a maximum  m one  study  and  a minimum  in 

^ the  observed  changes  in  water  level  during  the  melting  ibe  other  could  be  the  result  of  the  scaling  relationships 

period  indicate  that  this  was  a distinct  possibility  (sec  employed  to  relate  the  models.  Tlie  sc.iling  procedure 

i Fig.  13).  A reduction  in  the  heal  load  would  naturally  used  postulates  a proportional  relationship  between  the 
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tlow  rates  nl  the  Inso  iiiikIcIs,  lor  the  Model  I and  II 
tests  these  scaling  lelationships  predict  that  the  behavior 
ot  the  models  should  he  similar  it  the  Ilow  rates  are 
ssithin  a lalio  o(  2.5  to  I . Thus  the  .VO-gpm  lest  of 
Model  I should,  under  a scaled  heal  load,  he  the  same 
as  the  7.5-gpm  Model  II  lest  However,  the  test  results 
do  not  indicate  this:  instead  they  indicate  that  a non- 
linear relationship  may  exist  between  the  two  models. 

It  we  assume  that  a Ilow  rale  range  of  .VO  to  4.0 
gpni  in  Model  I actually  corresponds  to  the  4 75-  to 
7.5-gpm  /one  of  Model  II  (rather  than  the  predicted 
/one  of  7.5  to  10.0  gpm ) a similarity  in  results  emerges. 
Under  such  a relationship  both  models  would  then  ex- 
hibit an  average  outlet  water  temperature  which  in- 
creases to  a relative  maximum  at  low  Ilow  rales,  de- 
creases to  a relative  minimum  at  inlerinediatc  rales 
and  again  increases  pnor  to  levehng-off  at  even  higher 
Mow  rates  thereafter.  With  outlet  water  temperatures 
dependent  upon  both  Ilow  rate  and  density  inlluences, 
a linear  scaling  relationship  is  not  very  probable  and 
as  such  can  only  be  used  as  a guideline.  A method  of 
predicting  possible  nonlinearitics  in  the  Mow  rate  scal- 
ing is  discussed  in  the  section  which  presents  the  ap- 
proximation of  the  heat  transfer  crxrlTieient.  It  should 
be  emphasi/.ed  that  the  curves  in  Figure  1 1 are  not  to 
he  construed  as  exact  but  only  as  indicative  of  gen- 
eral vanations  in  temperature  between^different  flow 
rates. 

Although  little  information  is  available  concerning 
the  heat  and  mass  transfer  at  an  ice/water  interface 
under  conditions  exactly  comparable  to  those  of  the 


heal  sink,  the  woi  k ol  \'.imci  coikciiiiii'j  m.iiui.iI 
convective  melting  along  a vcilical  icc  ^ilccl  doc-' 
provide  some  insight  as  to  why  a imimmmi  icmpci.i- 
lure  miglil  occui  within  this  langc  ol  Ilow  laics 
Vanicr  noted  lhal  ihc  net  diicclion  ol  waici  iiiosc 
nicnl  was  upward  for  ceilam  langes  ol  bulk  w.iici 
leni|)eralurc.  a direct  lesull  ol  the  density  imcisioii 
of  water.  In  addition  he  observed  ihc  exisiciisc  ol 
/ones  where  the  net  motion  was  dowiiwaid.  ,is  well 
as  /ones  of  dual  Ilow  m which  walei  Mowed  m both 
directions.  Tlie  range  ol  these  dual  /ones  lends  to 
be  ralhei  narrow,  as  is  illustiated  m I iguie  I 2 which 
IS  reproduced  from  Ins  paper 

In  the  ice-watcr  heal  sink  the  gioss  Mow  iliieclion 
IS  always  downward;  lhal  is,  he.itcd  walei  is  added 
at  the  top  and  the  same  (|iiantily  is  removed  lioni 
the  bottom.  As  discussed  pieviously.  at  exliemely 
low  Ilow  rales  the  water  appears  to  move  dowiiwani 
in  a hori/onlal  layer  and  at  the  highest  Mow  laies  it 
moves  as  conventional  Mow  ihioiigh  an  aniuilai  sp.ice 
Tlicse  observations  indicate  lhal  a iraiisilion  occuis 
as  the  Mow  rate  increases.  The  density  nillneiices 
which  produced  stratirication  at  low  Ilow  rates  could, 
at  some  intermediate  rale,  provide  sufficient  resis- 
tance to  Ihc  annular  flow  priKCSS  to  generate  an  un- 
stable condition  where  the  motion  of  water  alter- 
nately traveled  upward  and  downward,  dependent 
upon  fluctuations  in  the  bulk  water  temperature. 

Such  a phenomenon  would  result  in  fluctuation  in 
the  outlet  water  temperature  similar  to  that  observed 
in  test  4,  namely  an  increase  in  temperature  up  to  the 
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mjviinimi  iloiisity.  tollowctl  h>  ;i  siiiKicn  tlccrcast.  iii 
iciiipcraliirc  as  Ilic  llnw  reversed. 

Xpproximalion  of  the  rale  of  melting 

rile  obsei'ved  change  in  water  level  with  lime  during 
eacli  expeinneni  tan  easily  be  cnnveiletl  lo  a gross  rate 
111  inelling.  Since  llie  lolal  mass  of  llie  sysiem  re- 
mains conslant.  llien 

whore  A/^  is  lire  mass  ol  ice  m sink  (Ibin).  and  Af^  is 
lire  mass  ol  walei  in  sink  (Ibm).  Tims 

,M  Mf,  JAf^ 

-Jl  = 0 = — If  — 

(//  (Jl  (ll 

Ml 

<//  til 

Similarly,  if  llie  ice  remains  submerged,  then  ilic  lolal 
volume  f,  IS  given  by 

where  k',  is  the  ice  volume  (ft’  ) and  is  the  water 
volume  (ft’ ).  Til  us, 

dt  dt  dt 
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•Since 

where  p,  is  the  density  ol  ice  and 
= Pw  > « 

where  is  the  density  ol  waiei.  ihen  lor  conslani 
densities 

‘Hj  =Ii^ 

dl  P| 

and 

dl  p^  dt 

Upon  subslituiion  tins  yields 

dt  Pj  dl  p^  dl 

Using  the  above  relationship  between  the  rales  of 
change  in  ice  and  water,  we  have 
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wlis'ls'  -1//  IS  cliiiiiiis'  III  ss.ils'i  Is'ss'l  sliiiiiii;  ills'  Unis'  in- 
sls'llls'lll  A/  Ills'ls'll'Is'.  ills'  Mis'  III  IIK'lllll^  S .III  Iv  .ip 
pioMlii.ils'sl  h\ 


Ills'  sll.lli;.'s's  III  ss.ils'i  Is'Ss'l  Ini  s'.lsli  Is'sl  ills'  slinssii 
III  Xpps'llilis  \ llUllis's  \l  \s  ll  sliniilsl  lis'  mils'll 
lli.il  ills'  ssuls'i  Is'scI  sls'sic.iss's  sliiiiii>;  iiu'lliiij;  muss'  is's' 
lias  .1  Invsci  sls'iisiu  lliaii  ssals'i  Oiks'  Ills'  is's'  lias  soiii- 
pls'ls'ls  iiis'lls'sl.  Ills'  vsals'i  Is'Vs'l  nisis'ass's  as  a Is'siill  nl  ills' 
sls's is'.iss'  III  ssals'i  Is'iisiis  vsilli  iiis'is'asiiij:  Is'iiips'ialiiis' 

III  Ills'  ahnvs'  sls'iisallnil  II  is  assiiiiis'sl  llial  ills'  sis'll- 
sils  nl  ills'  ssals'i  Is'ili.iiiis  sciisl.iiil , lliis  is  mil  pis'- 
siss'ls  Inis'  hs'saiiss'.  is  iiis'lliiij!  pinss's'sis.  ihs'is'  is  an 
ilKis'ass'  III  ills'  bulk  ssals'i  Is'ilips'ialiiis'  llnsssss'l. 

Illis  Is'ilips'ialiiis'  IlKis'ass'  Is  small  ami  ills'  assiiiiiplinii 
mlinslussrs  iiiilv  a small  s'lini  Ills'  bulk  ssals'i  Is'iii- 
ps'ialiiis'  /,,  IS  appinximals'il  as  ills'  nums'iis.il  avs'i.ijis' 

1)1  Ills'  iiils'l  ami  niillcl  ssals'i  Is'intH'ialuis's.  Ilial  is 

V 

sslis'is'  II  IS  ills'  Inlal  mimbs'i  nl  ts'Ui|K'ialuis'  rs'aslmns 
laks'ii  linm  tbs'  sl.iil  nl  ills'  Is'sl  iiiilil  sniiipls'linn  nl 
lUs'llmp. 

Allliniigli  s'i|  7 is'piess'iils  a lallis'i  siiiils'  appinvima- 
linii  ol  llic  bulk  ssals'r  is'iiips'iatiirs'.  bugs'  sluuigs's  in 
Ills'  bulk  water  Icmperalurc  rcsull  in  niily  small  eirnrs 
in  llie  rale  of  ice  melling.  l or  example,  if  llie  bulk 
water  temperature  were  calculatesl  to  be  40“F  but  wa.s 
actually  the  density  factor  would  cliange  be  a 
mere  0.6%.  In  Ibis  example  the  ice  density  was  taken 
as  56.7  Ibm/fi’,  rather  Ilian  the  standard  57.2  Ibm/fl* 


foi  sleai  ICC.  bs'caiisc  ibc  air  bubbliiii;  sysicm  caiivs'd 
s'liliamms'iil  nl  a large  iiiimbs  i nl  .m  bubbles  m llie  ice 
and  ibus  lediiced  Us  slensils  ( I leuic  '>  slinsss  an  bub- 
bles released  diiniiu  melliiig. ) Ifs  snmpaiiiie  ibe  milial 
ice  vobiine  prior  lo  ans  melline  in  ibe  lolal  sbaiiee  iii 
volume  loi  ibe  sink,  il  ssas  deleimiiied  ib.il  Ibe  slensils 
ol  Ills' ice  vans'll  Iroiii  .“in.f)  Ibiii  ll ' lo  .sb.')  Ibin/ll  ’ 
loi  Ills'  six  Is'Sis.  ills'  overall  aseiage  value  being 
Ibiii  I I ' 

Using  a lO-boni  mieival  in  ibe  piecediii'j  appioxima- 
iioii  s'liualioii.  Ills'  nils'  nl  meinn;'  ss.is  deleimmed  Ini 
each  lesl.  Ibese  lesiills  aie  plolled  in  Appemli.x  l ig- 
uies  AI-A5  and  labulaled  in  Table  l\'  'Mso  incinded 
III  lable  IV  IS  a ginss  inellme  nils'  ilelerimiied  bs  divid- 
nig  ills'  Inlal  iiiilial  Ke  mass  bs  ibe  lime  leqnireil  in 
Ills'll  all  Ills'  ice.  .\  cnmpaiisoii  nl  ibe  gmss  melling 
i.ile  lo  ills'  aveiage  value  ol  ibe  .ip[iioxuiials'il  meliiiig 
nils'  sliosss  a bigli  slegiee  ol  coirelalinii. 

A comparison  of  ibe  luellmg  rales  lor  ills'  dirieieiil 
lesls.  as  111  l igiiie  l .n  is  complicalod  bs  ibe  Iasi  Ibal 
Ills'  applied  Ileal  reiection  lale  vailed  loi  each  lesl  (sec 
Table  III  Tni  example,  in  conveiling  ibe  nils'  ol  niell- 
nig  iiilo  a rale  of  beal  exebange  (i.e.  each  pound  iiiellesl 
n's|uiied  144  Hiii),  ibe  beal  leieclion  lales  for  ibe  liisl 
l()-boui  inlersals  ranged  fiom  .U  .5.)(i  Ulu/bi  lo4I.IM 
Blii/br  Tliese  values  represenl  onis  llie  lalciil  beal  por- 
lion  ol  ills'  lolal  beat  added  lo  llie  sink; some  heal  isslorcd 
by  sensible  bealliig  ol  Ibe  svalei  imlially  iii  ibe  annulus 
and  ills'  iiielled  svalei  geneialed  ilnimg  llie  mleival. 

Tils'  perceiilage  of  sloied  sensible  beal  sanes  svilb  lime 
lioiii  only  aboiil  2'/r  during  ills'  early  slages  ol  heal  le- 
lesiion  up  lo  aboul  50'^  diirmg  Ibe  laler  stages  ol 
niellmg.  ll  sliould  also  be  noted  thal  tbs'  walei  lable 
elevations  in  Ibe  lank  were  measuied  wilbin  '/i„  in  of 
ills'  Inis'  value,  so  ibal  a maximum  eiioi  of  Vh  in.  oi 
appioximals'ly  I'/i  ssas  possible. 

In  reviewing  Ills'  data  preseiiled  in  Table  IV  (ex- 
cept for  ills’  results  of  lesl  4)  all  niellmg  rales  be 
ssitbui  a loleiance  band  of  approximalely  lO'v  of  llie 
aveiage  rates.  Alter  .^0  bours  of  operalion  lliere  is  an 
even  smaller  varialion  in  nicli  rale  among  llie  various 
lesls,  indicating  ibal  How  rale  effects  are  not  signili- 
cunl  in  the  litial  bours  of  Ills'  melting  process,  ll  is 
during  tins  |K'nod  that  the  remaining  ice  exists  m a 
small  mass  ai  the  lank  centerlme.  pioduciiiga  large 
annular  passage  for  water  How  and  lliiis  a subslanlial 
reduction  in  the  water  velocity.  Tlie  How  rale  inllu- 
eiicc,  manifested  by  Hie  ice  profiles  (Fig.  X),  exerts  a 
decreasing  iiiHuence  a.s  more  ice  is  iiielied,  until  a 
point  is  reaclied  when  all  sinks  begin  to  exhibit  llie 
same  basic  melting  rate. 
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Table  IV.  Melting  rates  - Model  II. 


Mvltim: 

Mrltiny 

rate  Jttnni' 

1»  V >•! 

<,rt)%\  inell 

Test 

f'lttw  rate 

rinu' 

rinfc 

tnft'tx'uf  ( 

f(  iniafmii 

appnn  rate 

(fiy'  rate* 

nti. 

fill) 

it  /O 

lit  .211 

2it-  Ui 

m All 

40- .50 

.50  6 0 

ilhm  hrt 

tihni  hri 

1 A 

4 7 \ 

63 

244 

242 

2 3s 

220 

|4K 

176 

220 

2 10 

IB 

A 

6S 

2K6 

264 

227 

20S 

141 

161 

222 

204 

2 

2 ^0 

S6 

2 U. 

22H 

22K 

206 

1 K4 

* • 

2 16 

2 14 

3 

lOOU 

5H 

2K6 

271 

242 

21  3 

176 

• » 

23H 

220 

4 

7 so 

67 

2 14 

20S 

IH3 

176 

164 

147 

142 

143 

S 

1 3 SO 

hO  h.l'^ 

244 

227 

227 

2 1 ^ 

14  1 

147 

204 

218 

•\vorago 

\alues^*^ 

261 

246 

232 

2 t 1 

1 MK 

161 

2 16 

2 1 1 

• c 

rnvs  melting  rale 

Initial 

ice  mass 

fol.il  ntcihng  linu' 


+ I’recisf  V aluc  uii  known  due  lo  equipmenC  failure,  me.in  v alue  of  61^  futufv  used  !•»  eaU  ulale  grttss  melt  in  r.ite 
••  Meilintr  completed  during  interval, 
tt  Average  excluiling  leM  4 


f'igurc  13.  .Mi’lliHf;  raws  fur  Model  // 


Tlic  variations  in  melting  rate  dunng  the  initial 
hours  probably  result  primarily  from  different  rates 
of  heat  refection  and  different  initial  ice  temperatures, 
rather  than  from  differences  in  flow  rates.  During 
this  penod  less  than  5%  of  the  rejected  heat  is  stored 
as  sensible  heat  , thus  the  rate  of  melting  is  almost 
identical  to  the  rate  of  heat  rejection.  For  example, 
the  differences  in  melting  rale  for  the  two  4.7.^-gpm 
tests  result  from  different  initial  ice  temperatures.  In 
test  I A the  initial  average  bulk  ice  temperature  was 
1 7. 6’ F,  while  in  test  I B it  was  3 1 .4° F.  As  indicated 
by  Figure  13  and  Table  IV,  the  test  with  the  lower  ice 
temperature  stored  a greater  amount  of  its  rejected 
heat  in  the  sensible  heating  of  the  ice.  Yet  both  tests 
show  a parallel  relationship  during  the  later  hours  of 


operation,  indicating  that  the  melting  laie  cveiituallv 
becomes  independent  ot  the  initial  conditions.  In  a 
similar  manner  the  variations  depicted  in  the  ol!  er 
tests  during  the  initial  hours  most  likely  result  liom 
different  heat  loads. 

Tlie  results  of  the  7.50-gpni  test  tend  to  ilhistiaie 
a possible  reduced  heat  rejection  rate  Although,  lor 
this  test,  the  average  heat  rejection  rates  shown  in 
Table  II  do  not  indicate  a low  heat  rejection  rate,  the 
very  low  rate  of  ice  melting  diK’s  imply  this  to  be  the 
case. 

Using  the  average  melting  rate  data  from  Table  IV, 
it  is  possible  to  show  the  partitioning  between  the 
amount  of  heat  rejected  to  the  ice  and  to  the  water 
for  each  lO-hr  interval.  Because  of  its  much  lower 
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moll  TJti's.  lest  4 (Mmlel  11-7.50  ppni)  will  bo  excliidcil 
Iruni  iho  avciage.  Diiriii};  llio  first  10  hours  the  average 
gross  rate  oT  melting  is  2()l  Ibm/lir.  With  the  rate  of 
melting  known,  the  rate  of  heat  storage  assoeiatcrl 
with  melting  can  be  determineil.  smee  144  Btn  is  rc- 
c|uired  to  melt  each  ponml  of  ice.  Using  both  the 
nominal  applied  heal  load  of  42,01. 7 Blii/hr  and  the 
average  rale  of  .17.700  Blu/hr  shown  in  Table  II.  the 
[X’rcentage  of  liilal  heal  involved  in  ice  niellmg  can  be 
denved.  Table  V lists  results  of  this  calculation. 

Tliese  results  indicate  that,  on  the  average  diiniig 
the  initial  lO-hr  period  (representing  about  KiTf  of  the 
total  melting  lime),  between  dO'f  and  of  the  re- 
lecled  heal  is  involved  m nicllmg  ice  During  the  last 
lO-hr  period,  an  average  of  between  55%aiul  00',7  of 
the  rejected  heat  Is  involved  with  ice  melting. 

If  How  rale  cflecis  exist,  they  apparently  do  not 
have  any  significant  intluencc  on  the  overall  rate  of 
niellmg.  Rather,  the  elfeci  ol  How  rale  is  manifested 
by  the  outlet  water  leni|X’ralure  dunng  melting,  a re- 
sult of  the  degree  of  overall  mixing  of  the  water  as  it 
passes  through  the  sink.  Tims,  for  low  How  rales  the 
slight  degree  of  mixing  ol  the  inlet  water  with  the 
melt  water  results  in  low  outlet  water  temperatures. 

As  flow  late  is  increased,  mixing  elfects  nicrea.se.  re- 
sulting in  higlier  outlet  water  temperatures. 

^.M  V.  Rates  ol  heal  rejei  led  to  the  ice  texclndinu  lest  4). 


nmv 

pentnJ 

jhri 

Rate  of 
meZ/iny  • 
ithm,hr} 

(HW  hr) 

(%) 

V,  V/' 
r.) 

0 10 

261 

37.6  13 

89.5 

99.5 

10  20 

246 

35.482 

84.5 

94.0 

20  30 

232 

33.379 

79,4 

88.3 

30  40 

21  1 

tO.442 

72.5 

80.6 

40  50 

1 MR 

27.072 

64.4 

71.6 

50  60 

161 

23.232 

55.3 

61.5 

* [>ala  from  Tablt*  IV. 

' 144  ./A/j. 

**  C^tiom  ‘*2.013  lliu/hr  (nominal  value). 
t+ 37.760  Hlu/hr  (Table  II). 

( OMPARISON  Of  EXPERIMENTAL  AND 
C OMPUTED  RE.SULTS 

As  miportani  as  the  determination  of  the  operating 
characteristics  of  the  heat  sink  may  be.  it  is  of  equal  im- 
portance to  determine  how  accurately  these  character- 
istics can  be  predicted.  A computer  program  developed 
dunng  the  Model  I sludy^  was  used  to  make  these 
predictions.  In  this  report  the  following  factors  are 
used  as  a basis  for  comparison:  mean  outlet  water 
temperature  during  the  melting  period,  overall  time 
required  to  complete  melting,  and  behavior  of  the 


outlet  water  temperature  after  all  the  ice  is  melted. 
I>ie  to  the  possible  variations  in  heat  rejection  rate 
for  the  different  tests,  two  computer  solutions  were 
developed  for  each  test.  One  used  a heat  rejection 
rale  e(|ual  to  the  value  determined  by  the  previously 
nieiilioricd  Prixedure  I involving  the  measurement 
of  the  coolant  water  How  rate  and  the  lenipeialure 
difference  mamlaincd  between  the  inlet  and  the  out- 
let. Tlie  other  solution  used  the  heal  rejection  rate 
lor  the  lower  of  the  values  determined  either  by  Pro- 
cedure 2 (the  heat  balance  technique)  or  Procedure  .1 
(the  slope  of  the  outlet  icmperalure-linie  curve.  This 
established  the  most  probable  band  of  heal  rejection 
rates  for  the  actual  lest. 

Wlieii  the  computer  solution  is  calculated  u-mg 
the  mniiniuni  heal  rejection  rale  as  determined  by 
Procedure  2 or  .1,  it  is  necessary  to  adjust  either  the 
input  coolant  water  flow  rale  B'or  the  temperature 
difference  across  the  sink  AT"  from  the  nominal  ex- 
[K’ritiienlal  values.  Tlie  How  rate  was  selected  for  ad- 
juslnicril  because  the  larger  number  of  temperature 
readings  taken  during  the  tests  led  to  a higher  degree 
of  confidence  in  these  results.  Tlie  adjustment  was 
made  to  the  How  rale  by  use  of  the  equation: 

,.|„w,ate  = -^2iSi-=— (gpm).  (8) 

fp^.AT  500. 1 CAT 

Table  VI  sunnnari/.es  the  heal  rejection  rates  used  for 
the  vanous  computer  runs.  Tlie  adjusted  How  rates 
are  listed  in  Table  Vlb. 


Table  VI.  Summary  «f  computer  predictions. 

Mean  temp 


Test 

no 

Heat  re/et  • 
titot  rate 
(fituthr) 

h'low  rate 

(gpm) 

during  melting 

inlet  outlet 

Melting 

time 

(hr) 

a.  Maximum  heat  rejection  rate 

) 

1 A 

41,305 

4.73 

43.6 

60.8 

61.5 

IH 

40.150 

4.73 

43.3 

60.0 

62.5 

2 

41,140 

2.50 

40.2 

72.7 

58.5 

3 

41,4)0 

10.00 

45.8 

53.9 

57.5 

4 

42.090 

7.50 

44.9 

55.2 

56.5 

5 

40,260 

13.50 

46.1 

52.0 

61.5 

h.  Minimum  heat  rejection  rate  (C?min 

) 

1 A 

38,600 

4.41* 

42.3 

59.6 

64.5 

1 1) 

35,050 

4.  1 2* 

41.8 

58.6 

70.5 

2 

40.720 

2.43* 

39.8 

72.3 

60.5 

3 

38.530 

9.30* 

45.0 

53.1 

6t.S 

4 

35.920 

6.40» 

43.3 

54.4 

67.5 

5 

34.790 

I2.66* 

45.3 

51.2 

64.5 

* Row  rate  adjusted  using  eq  R. 


In  Appendix  Figures  A I -AS  the  computer-pre- 
dicted  outlet  water  temperatures  have  been  plotted 
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with  the  comparable  experimental  test.  Inlet  tem- 
perature curves  for  the  computer  solutions  have  been 
eliminated  in  these  figures  since  they  dif  fer  from  the 
outlet  by  a constant,  the  AT" show'll  in  Table  III. 

Table  VII  and  F-igiire  I la  illustrate  the  variation 
between  the  computer  predictions  and  the  experi- 
mental observations,  (k-nerally.  the  computer  version 
predicts  a higlier  mean  outlet  water  temperature  dur- 
ing the  melting  priKcss  for  both  heat  rejection  rates. 

It  had  been  expected  that  the  use  of  the  two  heat  rc- 
lection  rates  miglit  result  in  the  predicted  tempera- 
tures bracketing  the  observed  values.  Althoiigli  the 
computer  program  did  not  predict  the  inllection  m out- 
let water  temperatures  which  iKCurred  in  the  neiglibor- 
hood  of  7.5  gpni.  it  otherwise  reasonably  replicated 
the  overall  increasing  trend  in  outlet  water  temperature 
wtth  increasing  coolant  water  flow  rate.  Tliis  indicates 
that  the  computer  model  simulates  the  gross  pr(x;ess 
well  but  predicts  a sliglitly  greater  amount  of  mixing 
of  the  inlet  water  and  the  water  within  the  sink. 

Table  Vli.  Average  outlet  water  temperature  during 
melting:  experimental  results  and  computer  predictions. 


Frediited  Predicted 


Test 

mt 

Measured 

Jn>m  Qmax  * 
(FI 

from  IJ„„,  • 

( 'F) 

1 A 

40.8 

43.6 

+ 6.7 

42.3 

+ 3.5 

1 1) 

41.2 

43.3 

+ 4.9 

41.8 

+ 1.4 

2 

38.5 

40.2 

+ 4.5 

39.8 

+ 3.5 

3 

43.4 

45.8 

+ 5.4 

45.0 

+ 3.4 

4 

40.3 

44.9 

+ 1 1.3 

43.3 

♦ 7.4 

5 

43.6 

46.1 

+ 5.8 

45.3 

+ 3.9 

Average  error 

♦ 6.4*5? 

*3.9*;? 

Tins  tendency  to  predict  a higlier  degree  of  mixing 
IS  well-illustrated  in  the  low  flow  rate  tests  after  all 
ice  has  been  melted.  IXrring  the  experimental  tests  it 
was  observed  that  the  predotninant  feature  of  the  post- 
melting penod  was  the  density  controlled  flow  of 
water  throiigli  the  sink.  Tliis  was  manifested  by  the 
motion  of  convective  water  cells  in  the  heat  sink  tank. 
As  the  test  flow  rates  were  increased  it  was  noted  that 
the  density  influences  had  a decreasing  influence  on  the 
post-melting  period  unli’  at  7.50  gpm  the  effects  were 
completely  absent. 

However,  for  the  lowest  flow  rate  the  computer  pre- 
diction of  the  outlet  water  temperature-time  history 
sliows  only  sliglit  fluctuations  of  the  type  associated 
with  density  influences.  Tliese  predicted  fluctuations 


are  more  representative  of  the  type  found  at  the  next 
higlier  flow  rate  |4.73-gpm)  expeiimeiilal  lest.  This 
leads  to  the  conclusion  that  the  flow  process  used  in 
the  computer  program  errs  sliglitly  in  predicting  the 
density-related  flow  effects  and  lends  to  predict  a 
greater  amount  ol  mixing  of  the  watei  within  the  sink 
than  actually  iKcurs.  Tins  limilallon  in  the  accuracy 
ol  the  prediction  docs  not  present  a serious  drawback, 
since  during  this  period  it  is  rather  easy  to  predict  the 
behavior  of  the  outlet  water  by  using  the  simple  rela- 
tionships developed  in  the  previous  section. 

lixeepi  for  the  deviation  cited  above,  the  computer 
program  does  simulate  the  behavior  of  the  sink  with  a 
liigli  degree  of  accuracy.  Tlie  errors  which  do  occur 
lend  to  be  conservative,  as  they  predict  higlier  than 
actual  outlet  water  temperatures.  Wlieii  the  predic- 
tion is  coupled  with  the  operation  of  a simulated  con- 
denser module,  a somewhat  lower  system  capacity 
will  be  predicted  than  Is  actually  available.  Assuming 
no  other  scaling  effects  are  encountered  with  much 
larger  heal  sinks,  it  appears  that  the  computer  pro- 
gram could  be  used  to  predict  the  outlet  temperature- 
time  curve  for  a prototype  sink.  Tlie  actual  outlet 
water  temperature  would  be  sliglitly  lower  than  pre- 
dicted. resulting  in  a higher  overall  system  efficiency 
and  subsequently  a sliglitly  increased  overall  capacity. 
Tlie  magnitude  of  this  nuHlest  increase  in  capacity 
would  actually  depend  upon  the  sensitivity  of  the 
plant's  efficiency  to  the  coolant  water  temperature. 
Tlius  the  error  in  the  prediction  of  outlet  water  temp- 
erature can  be  considered  to  be  a small  inherent  safely 
factor  in  the  selection  and  sizing  of  a prototype  sys- 
tem. 

A comparison  was  also  made  between  the  predicted 
lime  required  to  melt  all  the  ice  by  using  both  the 
higli  and  low  heat  rejection  rates  discussed  earlier  and 
the  lime  measured  in  the  lest  model.  Tliese  results  are 
presented  in  Table  VIII.  Tlie  higli  heal  rejection  rale 
results  in  a predicted  melting  lime  which  is  2.29f  less 
tlian  that  observed,  while  the  lower  rate  predicts  a 
melt  time  which  is  5.1%  gieatei  than  that  observed. 
Since  the  use  of  these  limiting  values  of  heal  rejection 
rate  lends  to  bracket  the  experimental  results,  it  is 
reasonable  to  conclude  that  the  program  is  predicting 
the  ice  melting  times  rather  accurately.  Tlie  predic- 
tion of  a higlier  outlet  temperature  than  observed  can 
be  attributed  to  inability  of  the  computer  program  to 
provide  for  the  effect  of  water  density  inversion  on 
tJie  flow  prixtess.  It  can  be  concluded  that  the 
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Table  VIII.  Comparison  of  melting  times  of 
experimental  and  computer  models. 


1 


reft 

n<>. 

MeasumI* 

(hr) 

Predicted 

(hr) 

Predicted 

f''"”  Qmm 

(hr) 

1 A 

63 

61.5 

-2.4 

64.5 

+ 2.4 

IB 

65 

62.5 

-3.9 

70.5 

+ 8.5 

2 

56 

58.5 

♦ 4.5 

60.5 

+ 8.0 

3 

58 

57.5 

-0.9 

61.5 

+ 6.0 

4 

67 

56.5  - 

10.5 

67.5 

+ 0.8 

5 

bo-sat 

61.5 

0.0 

64.5 

+4.9 

Average  error 

-2.2% 

+ S.I% 

• Correcled  lo  account  for  initial  icc  temperatures  below  32'*  F . 
t Precise  value  not  known  due  to  equipment  failure;  mean  value 
of  61. S hr  used  to  calculate  percentage  error. 


computer  program  tends  to  predict  a greater  amount 
of  mixing  of  inlet  with  outlet  water  than  will  actually 
(Kcur,  the  differences  most  likely  resulting  from  the 
density  inversion  of  water. 

COMPARATIVE  ANALYSIS  OF  THE  MODEL  SINKS 

Tile  first  experimental  investigations  of  tlie  ice- 
water  heat  sink  concept  employed  a test  apparatus 
4 ft  in  diameter  by  6 ft  high  (Model  I).  Tliese  initial 
studies  indicated  that  the  ice-water  heat  sink  concept 
was  both  feasible  and  practical  and  merited  further  in- 
vestigation. Ttie  results  of  these  prior  studies  are  pre- 
sented in  references  3 and  14.  Tlie  scaling  relationships 
developed  by  Brown  and  Quinn*  were  used  to  establish 
die  experimental  heat  rejection  and  coolant  flow  rates. 
Since  the  comparison  of  the  model  sinks  is  based  upon 
diese  scaling  relationships,  some  discussion  of  their  in- 
herent limitations  is  in  order. 

Based  mainly  on  geometric  considerations,  the  scal- 
ing relationships  were  developed  to  ensure  similarity  in 
the  operational  characteristics  of  different  size  sinks. 
Tins  would  Include  parameters  such  as  the  rate  at  which 
heal  could  be  rejected  lo  the  sink,  the  relationship  be- 
tween tbe  heal  rejection  rate,  coolant  water  flow  rale 
and  outlet  water  temperature,  and  useful  life  of  the 
sink.  Mathematical  difficulties  involved  in  modeling 
dynamic  effects,  such  as  the  density  inversion  of  water 
and  boundary  layer  flow,  precluded  their  inclusion  in 
the  relationships.  It  was  felt  that  these  influences 
would  be  of  secondary  importance  compared  to  the 


general  heal  transfer  priKess  and  that  simplified  scal- 
ing relationships  could  be  employed  to  relate  the  models 
with  a reasonable  degree  of  confidence. 

In  comparing  the  two  models  the  expectation  was 
not  to  find  a perfect  correspondence  in  test  results  hut 
rather  similarities  in  behavior,  such  as  the  melting  pat- 
tern, general  teniperalure-time  history  of  the  outlet 
water,  and  predictability  by  use  ol  the  computer  pro- 
gram. Minor  differences  in  the  performance  ot  compar- 
able tests  in  the  two  model  studies  would  not  be  con- 
sidered a serious  limitation,  as  variations  in  the  initial 
test  conditions  were  sufficiently  large  to  prevent  exact 
correlations. 

In  the  following  discussions  the  “Model  I tests”  refer 
solely  to  the  l6,80.S-Btu/hr  tests  conducted  using  the 
4-ft-diam  model,  and  the  “Model  II  tests”  refer  to  the 
42.013-Btu/hr  tests  conducted  using  the  6-ft-diam  ap- 
paratus. Althotigli  iiiaiiy  of  the  results  developed  dur- 
ing the  earlier  Model  1 studies  are  presented  in  tins 
report,  the  reader  should  consult  references  3 and  14 
if  additional  information  is  needed. 

Tlie  technique  for  measuring  melting  times  dif- 
fered in  the  two  test  series.  In  the  Model  I study,  the 
melting  time  was  determined  by  monitoring  thermo- 
couples and  by  direct  observation.  Tliis  technique 
dix;s  not  give  the  exact  time  at  which  melting  was 
completed  but  rather  it  shows  the  earliest  time  at 
wliich  ice  was  no  longer  visually  observed  in  the  tank. 
With  an  interval  of  about  4 hr  between  observations, 
it  is  possible  that  melting  may  have  been  completed 
2 to  3 hours  prior  to  the  observed  “melting  titne  ” 

In  the  Model  II  studies  melting  was  monitored  by 
measuring  changes  in  the  water  level  and  by  tempera- 
ture measurements.  Interpolation  between  points  on 
the  water  table  graph  (Fig.  A1-A5)  results  in  a finer 
measurement  of  the  exact  melting  time.  In  addition 
the  temperature  measurements  show  a high  degree  of 
correlation  with  this  method.  During  the  final  hour 
of  melting  a small  /.one  of  low  temperature  water 
exists  in  the  vicinity  of  the  remaining  ice;  this  zone  is 
located  along  the  vertical  centerline  of  the  tank  at  ap- 
proximately Vj  of  its  heiglit.  Once  the  ice  is  com- 
pletely melted,  the  temperature  of  this  zone  rapidly 
increases  lo  the  mean  temperature  of  the  water  sur- 
rounding the  zone.  For  the  higher  flow  rate  tests  the 
occurrence  of  this  sudden  temperature  change  coin- 
cides with  the  abrupt  change  in  the  behavior  or  the 
water  level  graphs.  In  the  lower  flow  rate  tests,  where 


I 
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density  stiatitkalion  effects  arc  present,  the  nsc  in 
temperature  of  this  /one  ikcuis  sliglitly  after  the  in- 
llection  in  the  graph  of  the  water  level,  the  difference 
normally  being  I to  2 hours.  It  is  felt  that  the  Model 

II  melting  times  arc  known  to  within  an  hour  of  the 
actual,  representing  a relative  error  of  one  hour  in 
sixty.  In  the  Model  I studies  melting  times  could 
vary  by  two  to  ihiec  hours,  for  a relative  error  of 
two  hours  111  forty  (with  the  observed  melt  time  al- 
ways longer  than  the  actual). 

Such  vanations  in  procedure  result  in  differences 
not  only  in  the  melting  time  parameter  but  also  in 
other  parameters  used  to  compare  the  models.  For 
example,  in  calculating  the  average  outlet  water 
temperature  during  melting,  it  is  averaged  up  to  and  in- 
cluding the  temperature  at  which  melting  is  completed. 
Since  in  the  final  hours  the  outlet  water  temperature 
sliows  a continuously  increasing  trend,  the  use  of  a 
melting  time  later  than  the  actual  will  increase  the  aver- 
age value;  this  amounts  to  about  4 to  5%  for  each  two 
hours  added  Similarly  any  parameter  dependent  upon 
the  melting  time  would  be  affected  by  this  difference 

III  procedure. 

To  show  the  relationships  between  the  tests  and  in 
the  two  model  sinks.  Table  IX  illustrates  the  flow  rate 
scaling  which  is  based  on  the  work  originally  reported 
in  reference  3 and  repeated  in  Ap[iendix  Bof  this  report. 


Table  IX.  Flow  rate  scaling  for  Models  I and  II. 

(Scaling  after  ref.  3.) 


MiHiel  // 

Test 

Flow  rate 

Test 

Flow  rate 

no 

(tfpml 

no. 

(gpm) 

l!>|.  IS] 

I.8Q 

1 A.  IB 

4.73 

2S 

1.00 

2 

2.50 

3S 

4.00 

3 

10.00 

4S,.4S2 

3.00 

4 

7.50 

5S 

0.63 

• 

1.55* 

• 

5.40* 

5 

13.50 

* No  tests  conducted  at  these  values. 


S 


The  predominant  feature  of  the  low  How  rate  tests 
was  the  development  of  a stratified  (low  process  involv- 
ing the  convective  motion  of  water  cells  through  the 
sink.  F-vident  in  both  models  this  indicates  a direct 
similarity  in  the  flow  pattern  between  the  1 .00-gpm 
Model  I test  and  the  2.50-gpm  Model  II  test.  The  inlet 
and  outlet  water  temperature-time  histories  for  both 


tests  are  sliown  m Figure  14.  A time  scaling  factor  of 
1.5  IS  used  to  provide  a comparable  time  basis,  m other 
words  the  actual  times  for  the  Model  I test  are  multi- 
plied by  1.5  (see  App  B for  denvatioii  of  the  scaling 
factor)  to  convert  them  to  Model  II  times.  Both  tests 
have  a relatively  constant  temperature  penod  which 
ends  soon  after  melting  is  completed  At  this  point 
there  is  a sudden  increase  m outlet  water  temperature 
to  the  temperature  of  the  inlet  water  a few  hours 
earlier,  this  is  felt  to  be  due  to  the  arnval  of  the  water 
cell  at  iiie  outlet.  Althougli  the  tests  were  ended  soon 
after  melting  was  completed,  it  is  believed  that  with  a 
constant  AT niamtained  between  inlet  and  outlet,  the 
inlet  water  begins  to  form  another  convective  water 
cell  at  the  top  of  the  sink  which  travels  downward  m 
a similar  iiiantier  to  produce  a second  constant  temp- 
erature penod,  followed  by  another  sharp  nse  m out- 
let water  temperature  If  this  were  not  tnie,  i.e.  if  the 
water  began  to  store  the  rejected  heat  by  mixing,  the 
outlet  water  temperature  curve  would  have  approached 
a constant  slope  line  similar  to  the  higher  flow  rate 
tests. 

Overall,  the  outlet  water  temperatures  for  the 
Model  I and  II  low  flow  rate  tests  show  a higli  de- 
gree of  correlation,  the  average  values  for  the  melting 
periods  diffenng  by  less  than  l .5°F.  If  a scaling  fac- 
tor of  1 .5  is  used  for  the  time  axis,  both  sinks  reach 
the  same  temperature-time  coordinate  at  approxi- 
mately the  same  relative  percentage  of  their  operat- 
ing time.  Ttiis  tends  to  confirm  the  general  applica- 
bility ol  the  Scaling  factors  for  low  flow  rates  for  the 
two  sinks  studied. 

Pie  next  higtier  flow  rate  tests,  those  of  1.89  gpm 
in  Model  I and  4.73  gpm  in  Model  II,  show  an  almos* 
exact  correlation  between  the  two  models  (Fig.  I 5). 
Pie  outlet  water  temperature  during  the  melting 
period  is  identical  for  the  two  models:  once  melting 
has  been  completed,  both  tests  experience  the  sliglit 
fluctuations  in  the  slope  of  the  outiet  curve  assrx'i- 
alcd  with  the  transition  from  a der'sity-infliienced 
flow  to  a mixed  flow. 

Puis,  in  the  lower  flow  rate  range,  the  scaling  re- 
lationships appear  to  provide  an  accurate  method  for 
relating  different  size  sinks;  however,  results  from  the 
intermediate  flow  rate  tests  indicate  that  these  rela- 
tionships may  not  be  used  for  all  flow  rate  conditions. 
Pie  results  of  the  intemiediate  flow  rate  tests  for 
3.00  gpm  in  Model  I and  7.50  gpm  in  Model  II  are 
sliown  in  Figure  16.  Fxcept  for  the  initial  hours  of 
operation,  distinct  differences  in  sink  behavior  iKcur. 
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Figure  15.  Comparison  of  4.  73-gpm  Model  II  test  to  1.89-gpm 
Model  I test. 
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Figure  16.  Comparison  of  7.50-gpm  Model  II  test  to  3.00-gpm 
Model  / test. 

With  two  tests  conducted  at  the  3.00-gpm  rate  dur-  gpm  to  10.0  gpm  determined  from  the  scaling  relation- 

ing  the  Model  1 study,  the  rather  sudden  rise  in  out-  ships  in  Appendix  B.  A more  detailed  discussion  of 

let  water  temperature  at  the  1 5th  hour  (Model  II  possible  nonlinearities  in  the  scaling  will  be  given  in 

time)  during  the  melting  period  was  substantiated.  the  section  concerning  the  approximation  of  the  heat 

As  noted  earlier,  the  possible  low  heat  load  during  transfer  coefficient. 

the  melting  periixi  in  Model  II  could  have  caused  Although  a comparable  high  flow  rate  test  was  not 

the  low  outlet  water  temperature  observed.  conducted  during  the  Model  I study,  the  13.50-gpm 

A comparison  between  the  highest  flow  rate  tests,  test  of  Model  II  indicates  an  interesting  trend  in  aver- 

4.0  gpm  in  Model  1 and  10.0  gpm  in  Model  II,  is  il-  age  outlet  water  temperature  during  the  melting  period, 

lustrated  in  Figure  17.  An  almost  constant  differ-  Between  the  interval  of  10  0 gpm  and  13.50  gpm  the 

ence  in  outlet  water  temperature  exists  between  average  outlet  water  temperature  becomes  essentially 

these  two  tests.  This  difference,  on  the  order  of  one  constant  at  a value  of  43.5°  (Fig.  I la).  Tliis  does 

to  two  degrees,  indicates  that  the  same  general  behav-  not  mean  that  during  the  melting  period  the  outlet 

ior  is  occurring  in  each,  but  that  a slight  difference  ex-  water  temperature  is  constant  (Fig.  1 7),  but  that  the 

ists  either  in  the  partitioning  of  the  rejected  heat  be-  average  value  over  the  entire  melting  interval  approaches 

tween  the  ice  and  water  or  in  the  overall  amount  of  a constant.  Should  this  trend  continue  for  even  higher 


inlet  water  mixing  with  the  water  in  the  sink.  This  is 
also  illustrated  in  Figure  1 1 which  compares  the  aver- 
age outlet  water  temperatures  during  melting.  In  the 
Model  I tests  there  is  a decrease  in  average  outlet 
water  temperatures  between  3.0  gpm  and  4.0  gpm, 
while  for  the  Model  II  tests  the  comparable  How  inter- 
val from  7.5  gpm  to  10.0  gpm  exhibits  an  increase.  In 
turn  this  indicates  that  the  scaling  relationship  between 
the  models  for  this  flow  rate  interval  is  not  correct  hut 
requires  additional  factors.  The  data  for  the  two 
models  would  exhibit  a higher  degree  of  correlation  if 
the  flow  rate  interval  of  3.0  to  4.0  gpm  in  Model  I 
were  scaled  to  correspond  to  an  interval  of  5. 5 gpm  to 
8.0  gpm  in  Model  II,  rather  than  the  interval  of  7.5 


fiows,  then  modeling  of  the  prototype  sink  and  con- 
denser unit  IS  simplified.* 

* The  similarity  between  the  outlet  temperature-time 
curves  for  the  10.00-  and  13.50-gpm  tests  is  depicted 
on  Figure  10.  The  somewhat  greater  ice  mass  in  the 
13.50-gpm  (Table  I)  is  responsible  for  the  difference 
in  melting  times.  Should  this  relation  hold  for  even 
higher  flow  rates  the  effect  of  flow  rate  on  the  rela- 
tionship between  outlet  temperature  and  time  could 
be  ignored  and  a "standard  "outlet  water  temperature- 
time curve  could  be  prepared  independent  of  flow 
rates  above  the  lO.OO-gpm  level,  or  the  appropriately 
scaled  value  for  another  size  heat  sink. 
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Figure  1 7.  Comparison  of  10. 00-gpm  Model  //  test  to  4. 00-gpm 
Model  / test. 


Comparison  of  Che  Model  1 experimental  results  to 
those  predicted  by  the  computer  program  is  dis- 
cussed in  reference  14.  In  general  the  same  type  of 
variations  observed  between  the  Model  II  experiment- 
al and  computer  results  also  occurred  in  the  Model  I 
study  (see  Fig.  I lb).  The  computer  program  tends 
to  predict  outlet  water  temperatures  slightly  higher 
than  those  observed  and  thus  represents  a conserva- 
tive prediction  for  both  Model  I and  Model  II  oper- 
ating characteristics. 


APPROXIMATION  OF  THE  HEAT 
TRANSFER  COEFFICIENT 

From  the  first  conceptual  studies  of  the  ice-water 
heat  sink  through  the  most  recent  experimental  test 
programs,  it  has  been  recognized  that  information 
relating  the  heat  transfer  coefficient  to  the  water 
temperature  and  flow  conditions  of  an  annular  flow 
sink  would  be  necessary  to  accurately  predict  the  be- 
havior of  a prototype  installation.  Mfith  little  infor- 
mation available  concerning  the  heat  and  mass  trans- 
fer at  an  ice/water  interface  under  conditions  similar 
to  those  in  the  heat  sink,  a series  of  progressively 
more  rigorous  methods  have  been  made  to  approxi- 
mate the  heat  transfer  coefficient  by  using  the  experi- 
mental data. 


The  first  method,  reported  in  reference  14,  used  a 
very  simplified  geometric  analysis  of  the  melting  pro- 
cess by  assuming  that  the  ice  cylinder  melted  uni- 
formly along  the  entire  vertical  surface  with  a con- 
stant rate  of  change  in  radius.  To  permit  a compari- 
son of  the  results  of  the  previously  utilized  approxi- 
mation with  the  refined  method  developed  in  this  re- 
port, a brief  discussion  of  the  earlier  techniques  is 
presented.  A rather  simple  relationship  was  devel- 
oped to  estimate  the  heat  transfer  coefficient  by  us- 
ing test  data  from  the  Model  I tests.  In  functional 
form  the  relationship  developed  (see  App.  C for  deri- 
vation) was 


^ _ Pj^  AR 
m m 


(9) 


where  h = average  heat  transfer  coefficient  (Btu/hr 
ft’  °F) 

AR  total  change  in  radius  of  tie  ice  cylinder 
(ft) 

= melting  time  (hr) 

= mean  temperature  difference  between  the 
ice  and  the  water  (°F) 

L = latent  heat  of  ice  (144  Btu/lbm). 


j 


It  must  be  cinphasi/cd  that,  in  this  approximation, 
the  rate  of  change  in  radius  was  assumed  to  be  con- 
stant and  equal  to  the  total  change  in  radius  (from 
the  initial  value  of  2 ft  to  ?ero)  divided  by  the  time 
required  to  melt  all  the  ice. 

Wlien  this  equation  was  employed  in  reference  14 
to  calculate  the  heat  transfer  coefficients,  the  melting 
time  used  in  the  relationship  was  not  the  observed, 
but  rather  an  effective,  melting  time.  It  was  reasoned 
that  sensible  heating  of  the  water,  rather  than  ice 
melting,  wxs  the  predominant  mode  of  heat  storage 
dunng  the  final  hours  of  melting  when  the  outlet 
water  temperature  began  to  increase  rapidly.  Thus 
it  was  assumed  that  ice  melting  had  essentially  term- 
inated at  the  lime  that  the  outlet  began  to  increase 
in  temperature  and  not  necessarily  when  ice  was  no 
longer  present  in  the  heat  sink.  For  the  low  coolant 
water  flow  rale  tests  the  effective  melting  time  agreed 
well  with  the  observed  melting  time,  but  in  the  higher 
flow  rate  tests  the  effective  melting  time  was  three  to 
five  hours  less  than  the  observed  values. 

Tlie  measurement  of  melting  rates  for  the  Model  II 
tests  casts  serious  reservations  on  the  use  of  this  effec- 
tive melting  time.  Unlike  the  assumed  behavior,  the 
actual  melting  rates  indicated  that  significant  amounts 
of  ice  are  still  being  melted  even  during  the  final  hours 
of  melting  (see  Table  II).  Tliese  results  indicated  that 
sensible  heating  does  not  become  the  predominant 
mode  of  heat  transfer  until  essentially  all  the  ice  is 
melted,  even  though  the  outlet  water  temperature 
has  begun  to  increase  before  that  time.  While  this 
might  seem  to  be  a contradiction  (i.e.  if  the  water 
temperature  is  increasing,  sensible  heating  must  occur), 
it  must  be  emphasized  that  the  behavior  of  the  outlet 
water  is  nut  synonymous  with  the  behavior  of  the  bulk 
of  the  water  in  the  test  tank.  Specifically,  sensible 
heating  is  manifested  by  an  increase  in  the  bulk  water 
temperature,  and  changes  in  the  bulk  water  tempera- 
ture are  not  necessarily  reflected  by  changes  in  the 
outlet  water  temperature.  This  behavior  is  most  evi- 
dent in  the  lowest  flow  rate  tests  where  the  outlet 
water  temperature  remains  constant  although  heat  is 
rejected  to  the  sink. 

The  effect  of  this  change  in  melting  time  of  the 
heat  transfer  coefficients  for  the  Model  I tests  is  illus- 
trated in  Figure  18.  Unlike  the  earlier  results  which 
indicate  that  as  the  flow  rate  increases  the  heat  trans- 
fer coefficient  increases,  the  revised  values  indicate 
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Figure  18.  Heat  transfer  coefficients 
for  Model  /,  16,805  Btujhr  tests. 


that  the  overall  heat  transfer  coefficient  tends  to  ap- 
proach a constant  value. 

The  revised  heat  transfer  coefficients  have  also 
been  tabulated  in  Table  X for  all  Model  I tests.  These 
values  should  be  used  instead  of  the  values  presented 
in  reference  14.  Also  reported  in  reference  14  was  a 
compilation  of  Nusselt  numbers  for  the  melting  pro- 
cess vs  the  test  flow  rate.  As  for  the  heat  transfer  co- 
efficient these  values  represented  overall  average  val- 
ues for  the  entire  melting  process.  In  calculating  the 
Nusselt  number  the  characteristic  dimension  selected 
was  the  diameter;  in  the  calculations  of  reference  14 
the  numerical  value  assigned  to  the  diameter  was  the 
initial  value  of  4 ft.  As  will  be  shown  later  in  this  sec- 
tion other  selections  for  the  characteristic  dimension 
can  be  made  which  result  in  a better  correlation  of  the 
data. 

Due  to  the  complex  melting  and  flow  processes  in 
the  heat  sink,  certain  limitations  or  restrictions  are  in- 
herent in  any  set  of  assumptions  employed  to  deter- 
mine heat  transfer  coefficients.  In  this  first  approxi- 
mation technique  the  major  limitation  or  discrepancy 
OCCUR  as  a result  of  the  assumption  that  the  rate  of 
change  in  radius  remains  constant  as  the  ice  cylinder 
melts.  To  illustrate  the  influence  of  this  assumption 
on  the  melting  process,  it  is  fairly  simple  to  develop 
an  expression  relating  the  change  in  radius  to  time. 

In  this  example,  nominal  values  from  the  Model  11 
tests  will  be  used,  since  the  rates  of  melting  for  these 
tests  were  measured  and  can  be  compared. 
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Table  X.  Heat  transfer  coefflcients  for  Model  I. 


Test 

no. 

Flow  rate 

(ft  pm} 

Time  to 

meit  all 

the  ice 
/hr/ 

Hulk 

wa  ter 

temp 

n-) 

Mean  temp 
Jiff  between 

ice  and  water 

Heat 

refection 

rate 

(Btufhr) 

Heat 

transfer 
coefficient 
(Btu/fH  hr^'h'} 

5S 

0.62 

36 

69.5 

37.5 

20,032 

1 2.00 

7S 

0.66 

S6 

54.9 

22.9 

11,272 

12.71 

2S 

1.00 

43 

53.7 

21.7 

15.200 

17.48 

6S 

1.89 

78 

45.0 

13.0 

8,425 

16.13 

8S 

1.89 

35 

53.4 

21.4 

29,520 

21.79 

IS, 

1.89 

41 

48.0 

16.0 

16.720 

24.83 

iSj 

1.89 

40 

49.0 

17.0 

18,320 

23.93 

4S, 

3.00 

33 

48.6 

16.6 

16,848 

30.18 

4S2 

3.00 

36 

47.9 

15.9 

16,776 

28.97 

3S 

4.00 

40 

45.4 

13.4 

16,350 

30.47 

If  it  is  assumed  that  the  change  in  the  radius  of  the 
ice  cylinder  is  uniform  over  the  length  and  occurs  at 
a constant  rate,  then 

constant  = C, . 

dt  ' 

Tlie  functional  relationship  between  the  radius  and 
time  in  operation  can  be  found  from  direct  integra- 
tion. Tlius 

With  the  boundary  conditions  that  at  time  zero 
IS  equal  to  the  initial  value  Rg,  and  at  the  melting 
lime  i„,R(t)  is  zero,  the  constants  C,  and  Cj  can  be 
evaluated.  This  results  in  a solution  of  the  form 

^(0  = /?o(' --)•  (10) 

Since  it  is  also  assumed  in  the  approximation  that 
melting  occurs  only  radially  along  the  cylindrical  sur- 
face and  not  in  the  vertical  direction,  the  volume  of 
ice  melted  at  any  time  l can  be  found  from 

(ft’)  (II) 

wliere  Hg  is  the  initial  height  of  the  ice  cylinder  (ft). 

For  the  Model  II  tests  the  average  values  of  the  con- 
stants appearing  in  eq  10  and  1 1 are 


Rg  = 3.0  ft 

Hg  = 8.6  ft 

61.76  hr  (Table  IV). 

Thus 

«(/)=  3.0(1 -0.0160  (10a) 

and 

F„(/)=  243-27|/?(Ol’.  (11a) 

Tlierefore,  after  the  first  10  hr  of  an  average  Model 
II  test  the  radius  would  be  approximated  as 

/?(I0)=  2.51  ft 

and  the  volume  melted  would  be 
F„(10)=27.01(9)-(2.51)’l 
F„(10)=  72.9  ft’. 

To  estimate  the  rate  of  melting  the  rate  of  change  of 
volume  can  be  approximated  by 

^ = 72.9  ft’ 

dt  ~ A/  10  hr 

Thus 
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= 7.29  ftVtu 


atui  llieretore.  the  rate  ol  rnelting  WDuld  be 

‘AW.  _ 
dt  Jt 

wtiere  p,  is  56.7  Ibin/t't’. 

Using  these  values  lot  the  first  10  hours,  the  ap- 
ptoxinmled  rale  of  melting  would  be 


— !-  = (56.7X7.29)=  415  Ibm/fu 

dt 


Tills  approximate  rate  of  melting  is  significantly 
greater  than  the  average  rale  of  261  Ibtn/hr  observed 
for  the  first  10-hr  period  (Table  IV);  in  fact,  it  ex- 
ceeds the  observed  rale  by  a factor  of  1 .6.  The  dis- 
crepancy, as  slated  previously,  results  from  the 
simple  a.ssumplion  that  the  rate  of  change  of  the 
radius  is  constant. 

While  this  method  of  approximation  was  accep- 
table when  the  rales  of  melting  were  unknown,  an 
improved  relationship  can  be  developed  to  estimate 
the  heat  transfer  coefficient  now  that  the  melting 
rates  are  available.  Starting  only  with  the  assump- 
tions that  during  the  melting  stage  the  ice  remains  in 
a cylindrical  form  and  the  melting  occurs  only  in  the 
radial  direction,  the  rate  of  change  of  the  radius  as  a 
function  of  the  melting  rate  can  be  determined. 

For  a cylindrical  geometry,  the  volume  of  the  ice 
is  given  by 

y,U)=nHJR{i)]\ 

From  this  the  volumetric  rate  of  change  of  the  ice  is 


Therefore,  the  rate  of  change  of  volume  is 


dV  \u 

— ^ 2nH„R{i)  — 
dt  ^ Ar 


Since  this  equation  will  be  employed  to  approxi- 
mate the  rate  of  change  dunng  a specific  time  interval, 
f,  lo/j,a  mean  value  for /?(r)  will  be  employed.  Tlius, 
if  A/?  IS  the  decrease  in  the  radius  during  the  tune  in- 
terval from  /,  to  fj,  then  the  mean  value  of  the  radius 
dunng  this  interval  is  given  by 


RU)  = RU,}- 


wtiich  upon  substitution  into  the  relationship  for  vol- 
ume yields 


Ao[a?0,)-^]- 


^=2tr//o/?(r)^^. 


This  expression  is  quadratic  in  terms  of  A/?.  Solv- 
ing for  A/?  in  terms  of  the  other  variables  yields 

/ ^ A,  tAK 

A/?  =/?(/,  )± 

tt/Zq  dt 

Since  R{t)  must  decrease  in  any  interval  from  f , to 
and  decrease  from  the  initial  value  R^  to  zero  dur- 
ing the  entire  melting  interval,  the  negative  root  must 
be  used.  Additionally,  since 


dt  Pj  dt 


As  in  the  previous  approximation  technique  it  will 
be  assumed  that  the  rate  of  change  in  the  radius  can 
be  estimated  by 


„ / , At 

AR=R{t,}-  R(t.f-— . 

' V ' Pi»r//o  dt 


dR  ^ AZ? 
dt  Ar 


While  presently  the  length  of  the  period  Ar  is  unspeci- 
fied, it  is  obvious  that  as  the  length  of  this  period  de- 
creases the  accuracy  of  the  approximation  increases. 
At  the  limit  as  A/  approaches  zero,  it  will  be  exact. 


From  this  relationship  it  is  possible  to  determine 
the  change  in  the  radius  of  the  ice  cylinder  during  any 
time  interval  from  /|  to  Ij,  provided  that  both  the  val- 
ue of  the  radius  at  the  beginning  of  the  interval  and  the 
rate  of  melting  during  the  interval  are  known. 

To  illustrate  the  improvement  in  the  estimation  of 
change  in  radius  by  using  eq  1 2,  Figure  19  has  been 
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I 


• from  eq  C2. 
t From  Table  VIII. 
••  From  eq  1 2. 
ft  From  eq  Cl. 
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tli-vcliipcd  will)  data  frimi  llie  two  4 Td-ppm  ifsis. 
tost.')  I A and  I B.  Also  iiKliidod  in  ihis  fijiiiic  are  la- 
dial  nieasurcnienls  taken  at  various  stages  ol  iiielling 
of  the  ite  cylinder  Itoin  the  two  I .H‘>-gpin  Model  I 
tests.  The  measurements,  which  weie  made  at  Vj  and 

of  the  tank  elevation  m the  Model  I tests,  have 
been  averaged  to  given  an  overall  cluutge  m radius.  A 
comparable  time  basis  between  the  Model  I and  Model 
II  tests  was  developed  by  using  the  scaling  factor  as 
discussed  m Appendix  B.  All  values  of  the  radius  have 
been  normali/ed  by  division  by  the  initial  radius  val- 
ues, .VO  ft  for  Model  II  and  2.0  ft  for  Model  I. 

A.S  illustrated  m Figure  1‘)  there  is  a higli  degree  of 
correlation  between  the  approximated  values  and  the 
measured  average  values.  In  general  the  correlation  is 
best  during  the  early  melting  stages;  the  largest  devia- 
tions occur  dunng  the  final  hours.  Tlic  decrease  in  ac- 
curacy of  the  approximated  solution  is  to  be  ex[>ecled. 
since  as  melting  proceeds,  the  discrepancy  increxses 
between  the  assumed  and  actual  melting  geometry . 
Althougli  errors  lend  to  be  greater  during  the  final 
hours,  this  revised  approximation  technique  represents 
a siibstaniiai  improvement  over  the  earlier  method. 

Ttiis  method  of  approximation  also  results  in  a dis- 
tinct improvement  in  the  estimation  of  the  heal  trans- 
fer coefficient.  Not  only  can  the  term  in  cq  d 

be  approximated  with  a higher  degree  of  accuracy,  it 
IS  possible  to  estimate  the  heal  transfer  coefficient 
dunng  different  time  intervals  within  a particular  test. 

As  with  the  earlier  approximation  there  are.  how- 
ever, definite  limitations  on  the  application  of  these 
equations.  Tlie  most  basic  limitation  results  from  vari- 
ations between  the  assumed  and  actual  melting  geo- 
metry. In  the  model  sink  the  rale  of  change  of  the 
radius  is  not  only  a function  of  time  but  also  of  the 
heiglii  of  the  ice  cylinder.  Additionally,  melting  oc- 
curs on  the  upper  and  lower  horizontal  surfaces,  liven 
if  these  other  effects  were  included,  some  allowance 
would  also  need  to  be  made  for  the  point  of  introduc- 
tion or  melt  scour  patterns  assiK'iated  with  the  inlet 
header  manifold  arrangement. 

Althougli  not  readily  apparent,  there  is  another  im- 
plied limitation  in  these  equations.  Wliilc  the  accuracy 
of  the  approximation  will  generally  increase  as  the 
length  of  the  lime  interval  A/  Is  decreased,  there  is  a 
lower  boundary  on  the  minimum  si/.e  of  this  interval. 
Since  the  equations  depend  upon  the  determination 
of  the  melting  rate  which  itself  is  calculated  by  detect- 
able changes  in  volume,  they  can  only  be  employed  for 


fairly  significant  time  intervals.  In  the  conversion  of 
ice  to  water  the  change  in  density,  and  thus  volume,  is 
only  approximately  W.  Tlicorctically.  it  would  be  pos- 
sible to  measure  this  change  for  as  little  as  10  lb  of  ice 
melted,  but  in  the  actual  test  apparatus,  a much  larger 
change  is  required  to  ensure  reasonable  measurements. 
Additionally  this  is  also  infiiienced  by  second  order 
effects,  such  as  changes  in  density  of  the  water,  expan- 
sion and  contraction  of  the  test  apparatus  and  evapora- 
tion of  water  from  the  sink. 

A.S  a result  of  these  limitations  the  approximation 
technique  hxs  been  applied  to  10-hour  intervals  cor- 
responding to  those  for  which  the  average  melting 
rates  were  determined.  In  using  these  10-hotir  inter- 
vals the  mean  temperature  difference  between  the  ice 
and  water  is  calculated  by  using  the  procedure 
described  in  Appendix  C as  applied  to  the  specific 
tunc  interval.  A sample  calculation  of  the  heat  trans- 
fer coefficients  for  test  I B is  illustrated  m 1 able  XI. 

Tlie  appioximated  heat  transfer  coefficients  for  all 
the  Model  I’  tests  appear  on  l iguies  20a  and  20b. 
wheic  the  variations  have  been  plotted  vs  the  cal- 
culated bulk  water  temperature.  These  results  have 
been  divided  into  two  tlow  rate-related  regimes:  a 
low  fiow  rate  /one  depicting  relationships  for  How 
rates  of  7.50-gpm  and  less,  and  a high  fiow  rate  /one 
depicting  the  10.00-  and  I.T50-gpm  test  results.  Also 
included  m Figure  20  ate  the  results  of  Tkachev'*  and 
Vainer"'  who  studied,  respectively,  turbulent  and 
natural  free  convective  melting  of  ice. 

Ttiere  is  an  excellent  correlation  between  the  ap- 
proximated heat  transfer  coefficients  and  the  results 
of  Vainer  for  the  test  flow  rates  of  4.73  and  7.50  gpm. 
Tlie  lack  of  correlation  for  the  2.50-gpm  tests  is  be- 
lieved due  to  two  influences.  For  this  test,  melting  ik- 
curs  in  a predominantly  vertical  rather  than  radial  man- 
ner. thereby  changing  the  surface  area  over  which  melt- 
ing iKcurs  from  that  associated  with  the  assumed  radial 
melting  configuration.  Tlie  other  effect  influencing 
this  variance  is  the  predominance  of  the  density-related 
flow  which  produces  a different  flow  pattern  from  that 
observed  in  the  4.73-  and  7.50-gpm  tests. 

As  flow  rates  were  increased  to  10.00  and  13.50  gpm 
a different  trend  became  apparent.  Unlike  the  low  flow 
rate  results  which  indicated  a gradual  increase  in  heat 
transfer  coefficient,  it  appears  that  for  these  two  tests 
the  heat  transfer  coefficient  becomes  approximately 
constant.  Additionally,  since  as  the  heat  test  prcKeeds, 
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the  bulk  water  temperature  generally  increases  (see 
Table  XI).  these  higlter  flow  rate  results  tend  to  pro- 
ceed from  a correlation  with  the  results  of  Tkachev'* 
to  a correlation  with  the  results  of  Vanier."*  Iliis  in 
turn  would  indicate  a possible  transition  in  How  be- 
havitw  from  “pseudo-turbulent"  free  convection  to 
“pseudonatural"  free  convection.  Ttie  use  of  the 
prefix  “pseudo”  should  be  noted,  since  in  all  the 
tests  there  is  an  imposed  overall  forced  flow  condi- 
liori.  that  is,  water  is  added  at  the  top  of  the  sink  at 
a specific  rate  and  removed  from  the  bottom  at  the 
same  rate.  Of  course,  this  does  not  actually  require 
tliat  the  water  “fiow"  through  the  sink,  since  in  die 
lowest  flow  rate  tests  it  was  observed  that  “flow”  was 
actually  the  motion  of  a hori/ontally  stratified  layer 
througli  the  test  lank. 

To  examine  possible  correlations  of  the  heal  trans- 
fer coefficient  and  coolant  water  fiow  rale,  the  ap- 
proximated heal  transfer  coefficients  have  been 
plotted  versus  a Reynolds  number  based  on  an  equiv- 
alent diameter  (see  App.  Cfin  Figure  21.  Tliey  have 
also  been  tabulated  in  Table  Xlla  for  test  fiow  rales 
of  4.7,^  gpm  and  less  and  in  Table  Xllh  for  flow  ratas 
of  7.50gpni  and  greater.  As  indicated  in  the  figure, 
two  distinct  trends  can  be  discerned.  For  the  test 
flows  of  2.50  and  4.73  gpm  there  appears  to  be  a sig- 
nificant change  in  heat  transfer  coefficient  with 
changes  in  Reynolds  number.  Unlike  the  behavior 
indicated  m Figure  20a,  the  results  of  the  2.50-gpm 
test  appear  to  correlate  with  the  others  on  a flow  rate 
basis.  Another  difference  with  the  trends  indicated 
in  Figure  20a  is  that  the  7.50-gpm  test  results  do  not 
indicate  any  correlation  with  the  results  of  the  4.73- 
gpm  tests;  rather,  they  indicate  a transition  between 
the  results  obtained  with  the  lower  and  higlier  flow 
rate  tests.  Again,  the  results  of  the  higliest  flow  rale 
tests  indicate  that  the  heat  transfer  coefficient  tends 
to  approach  a constant  value  over  the  range  of  Rey- 
nolds numbers  applicable  to  the  experimental  tests. 

Tile  trends  illustrated  in  Figures  20  and  21  tend  to 
indicate  that  some  form  of  correlation  exists  which 
can  be  used  to  relate  the  heat  transfer  coefficient  to 
the  applicable  lest  conditions.  The  most  general 
form  would  be  a relationship  involving  the  Nussell 
number  as  a function  of  other  dimensionless  para- 
meters. Before  discussing  these  relationships  some 
background  information  miglit  be  useful. 

In  general  two  basic  correlations  can  be  established, 
as  illustrated  by  Figures  20  and  21 . The  first  would  be 


Table  XII.  Heat  transfer  coefficient  for 
various  Reynolds  numbers. 

Test  Reynolds  number  Heat  transfer  coefficient  h 
no.  ([iti4/hr  ft^  ^ F) 


a.  For  2.50-  and  4.7.3-gpm  tests 


2 

44.0 

1 1.6 

2 

47.4 

1 I.S 

2 

51. 1 

1 2.4 

2 

56.1 

12.3 

2 

62.2 

13.2 

lA 

75.9 

16.8 

IB 

76.0 

19.9 

1 A 

81.8 

17.0 

IB 

82.4 

19.4 

1 A 

86.4 

19.0 

IB 

87.9 

19.7 

lA 

94.7 

20.5 

IB 

97.1 

20.2 

1 A 

106.0 

22.5 

IB 

109.5 

23.5 

1 A 

130.3 

29.1 

IB 

145.7 

34.2 

b.  For  7.50- 

, 10.00-  and  13.S0-gpm  tests 

4 

I 14.7 

18.5 

4 

120.4 

18.3 

4 

125.6 

18.8 

4 

133.9 

18.9 

4 

146.2 

19. 1 

3 

151.0 

27.1 

3 

167.5 

23.7 

4 

I 70. 1 

22.6 

3 

181.9 

23.8 

5 

198.9 

26.3 

3 

202.0 

23.4 

5 

2 17.2 

2 1.9 

3 

233.6 

22.3 

5 

2 34.5 

22.5 

5 

256.4 

23.2 

5 

289.1 

23.0 

5 

350.7 

19.1 

typical  of  a natural  convective  problem,  the  variations 
in  heat  transfer  coefficient  and  thus  Nusselt  number  re- 
sulting directly  from  temperature-density  induced  in- 
fluences. In  parametric  form  this  is  typically  expressed 
as 

Nu  = /,(Gr  Pr) 

where  Gr  and  Pr  represent  the  Grashof  and  Prandtl 
numbers,  respectively.  On  the  other  hand  the  results 
of  Figure  21  indicate  a different  form  of  correlation, 
specifically: 


29 


Roc  Reynolds  Number 


Figure  21.  Heat  transfer  coefficient  vs  Reynolds  number.  Model  11. 


Figure  22.  Variation  in  Nusselt  number  with  Reynolds  number,  Model  //. 


Nu  = /j(Rc  Pr) 

wlierc  Rc  is  the  Reynolds  number.  Tlie  questions 
vvtiich  arise  pertain  to  the  parametric  form  wliicli 
would  best  result  in  predicting  the  dominant  influ- 
ence and  the  variables  which  should  be  selected  as  the 
characteristic  dimensions  needed  to  nondimensional- 
i/e  the  parameters.  If  a relationship  of  the  (Gr  Pr) 
type  is  a.ssumed.  normal  practice  dictates  the  use  of 
the  length  of  the  ice  cylinder  as  the  characteristic 
dimension.  If  the  (Rc  Pr)  form  is  selected,  correla- 
tions arc  generally  based  on  the  use  of  some  type  of 
diameter:  however,  there  arc  instances  where  a length 
factor  has  been  included. 

In  the  selection  of  the  parameters  it  was  decided 
to  use  a relationship  of  the  ( Re  Pr)  type  with  a Rey- 
nolds number  based  on  an  equivalent  diameter  and  a 
Nus-sclt  number  based  on  the  length  of  the  ice  cylin- 
der. file  use  of  an  equivalent  diameter  for  correlat- 
ing Hows  111  annular  passages  is  fairly  common,''  how- 
ever. only  relationships  concerning  forced  convection 
could  be  lound  in  the  literature.  Tlie  use  of  length 
in  tfie  calculation  of  the  Nusselt  number  for  correla- 
tion with  Reynolds  number  based  on  a diameter  is 
somcwliat  of  a break  with  traditional  approaches  but 
was  selected  for  tw  > specific  reasons. 

Tlie  first  reason  is  a rather  practical  one.  If  the 
equivalent  diameter  is  selected  as  the  characteris- 
tic dimension  for  the  Nusselt  number,  it  produces  ex- 
tremely large  vanations  in  the  calculated  values.  For 
example,  when  an  equivalent  diameter  is  used  in  test 
lA  the  average  Nusselt  number  for  the  first  10-hour 
period  would  be  only  14.2.  On  the  other  hand  for  the 
last  10-hour  period  it  would  be  384.9.  Tlius  for  a 
739f  increase  in  the  heat  transfer  coefficient,  from 
16.8  to  29.1  (see  Table  Xlla),  the  Nusselt  number 
would  increase  by  approximately  2710%.  Tliis  does 
not  appear  to  be  a very  reasonable  method  on  which 
to  base  a correlation.  If  instead,  the  length  of  the  ice 
cylinder  is  employed,  the  73%  increase  in  heat  trans- 
fer coefficient  results  in  the  same  proportional  in- 
crease in  Nusselt  number.  Mere  the  characteristic 
length  is  assumed  to  remain  constant  at  the  initial 
length,  which  is  consistent  with  the  assumptions  used 
in  the  denvation  of  the  approximation  equations  for 
die  heat  transfer  coefficient. 

The  second,  and  more  fundamental,  reason  is  illus- 
trated by  examining  the  relationship  between  free  and 
forced  convection.  Free  convection  can  be  considered 


to  be  the  limiting  case  of  forced  convection  wlien  in- 
fluences of  the  Reynolds  number  type  can  be  neglected. 
As  indicated  by  Table  XII  and  Figure  2 1 , the  Reynolds 
numbers  for  the  Model  II  tests  are  indeed  small  and 
tend  to  approach  the  /one  wherein  they  would  have  no 
effect  at  all.  As  stated  previously  for  the  lowest  flow 
rate  tests,  density  (i.e,  free  convective)  related  influ- 
ences predominate  even  under  the  imposed  force  flow 
conditions.  Ttius  over  the  range  of  flow  rates  tested, 
the  Model  II  heat  sink  apparently  is  influenced  by  both 
density  and  flow  velocity  parameters.  Since  the  nor- 
mal correlation  of  free  convection  of  vertical  cylinders 
employs  the  characteristic  length  of  the  cylinder  as  the 
cliaractenstic  dimension,  it  was  felt  appropriate  to 
also  use  this  dimension  in  calculating  Nusselt  numbers. 

Tlie  results  of  this  analysis  are  presented  in  Figure 
22.  As  in  Figure  21  there  are  two  distinct  /.ones,  the 
first  applicable  to  flow  rates  of  4.73  gpm  and  less  and 
the  second  to  flows  of  7.50  gpm  and  greater.  Again 
there  is  some  overlap  between  the  two  /ones,  indicat- 
ing some  transitional  effects.  Tliis  correlates  with  the 
observed  transition  in  the  behavior  of  the  outlet  water 
after  melting  that  (Kcurred  within  this  same  range  of 
flow  rates.  Using  a least  squares  fit  for  the  two  differ- 
ent sections,  the  following  relationships  have  been 
developed: 

Nuj  = 0.881(Re„pPr)“''”  (13a) 

for 

450  < Repp  Pr  < 1400 
and 

Nug  = 39.72(RtppPr)‘’^''*  (13b) 

for 

1400  < Repp  Pr  < 3500. 

In  the  section  of  this  report  which  compared  the  re- 
sults of  the  two  different  models,  there  appeared  to  be 
a deviation  between  the  linear  (low  rate  scaling  relation 
ships  postulated  in  Appendix  B and  the  observed  aver- 
age outlet  water  temperatures,  as  illustrated  in  Figure 
II.  It  was  noted  that  a better  correlation  between  the 
outlet  water  temperatures  for  the  two  models  would 
occur  if  the  flow  rate  /one  from  3.0  to  4.0  gpm  in 
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Model  I were  shifted  to  correspond  to  tlte  Model  II 
range  from  4.75  to  7.5  gpni  rather  than  the  linearly 
scaled  range  of  7.5  to  10.0  gpni.  By  assuming  vari- 
ous ice  diameters  for  these  Model  I and  Model  II 
tests  and  comparing  Reynolds  numbers,  it  is  possible 
to  explain  the  probable  reason  for  this  nonlinearity  in 
scaling.  As  the  diameter  of  the  ice  cylinder  decreased 
from  .^.5  ft  to  0.5  ft  in  the  3.00-gpm  Model  I test 
(nominal  tank  diameter  was  4.0  ft),  the  Reynolds 
number  (based  on  an  equivalent  diameter)  varied  from 
77.5  to  I 2*^.  For  the  4.00-gpm  Model  I test  the  same 
changes  in  diameter  resulted  in  a variation  in  Reynolds 
number  between  d8  and  164.  In  the  Model  II  test 
with  a flow  rate  of  4.73  gpm  the  variation  in  Reynolds 
number  is  between  76  and  130.  This  would  indicate 
that  on  a flow  rate  basis  a more  comparable  flow  re- 
lationship exists  between  the  3.00-gpm  Model  I and 
4.73-gpm  Model  II  tests  rather  than  the  linearly 
Scaled  3.00-gpni  Model  I and  7.50-gpm  Model  II  ex- 
periments. Similarly  the  variation  in  Reynolds  num- 
ber from  d8  to  164  for  the  4.00-gpm  Model  I test 
closely  corresponds  to  the  7.50-gpm  Model  test  which 
had  a Reynolds  nvimhcr  varvation  from  f f 5 to  170. 
Althougli  this  analysis  indicates  the  limitations  of  the 
linear  scaling  relationships,  they  are  extremely  useful 
tor  quickly  estimating  heat  loads  and  flow  rates.  How- 
ever. once  these  initial  parameters  have  been  estimated, 
an  analysis  of  the  Reynolds  number  type  would  be 
more  appropnate  for  establishing  practical  operational 
parameters. 

Of  course  equations  13a  and  13b  cannot  represent 
a complete  characterization  of  the  melting  process.  As 
stated  earlier  there  are  density-related  effects  which 
undoubtedly  should  be  included,  possibly  by  means  of 
a Orashof  number  term.  Without  considerable  data 
available  at  a constant  Reynolds  number,  it  is  impos- 
sible to  define  this  functional  relationship, however, 
for  the  lower  flow  rate  zone  (i.e.  Re[,|;  Pr  less  than 
1400),  it  appeared  that  increases  in  the  Grashof  num- 
ber resulted  in  small  decreases  in  the  heat  transfer  co- 
efficient. Tliis  is  reasonable  if  the  density-related  ef- 
fects (upward)  are  opposed  to  the  normal  flow  direc- 
tion (downward).  It  is  important  to  recognize  that 
while  the  relationships  developed  in  fliis  report  used 
certain  parameters  and  characteristic  dimensions, 
other  satisfactory  correlations  are  quite  possible.  It 
is  felt  that  other  variables  such  as  the  aspect  ratio  of 
the  ice,  i.e.  the  length  to  diameter  ratio,  could  also  be 
important.  Due  to  the  limited  scope  of  this  study  and 


the  generalized  nature  of  the  calculations,  a more  de- 
tailed analysis  is  not  possible. 

CONCLUSIONS  AND  RECOMMENDATIONS 

A major  concern  recognized  during  previous  ex- 
perimental studies  with  ice-water  heat  sinks  was  the 
potential  problem  in  scaling  results  from  the  small 
model  to  a prototype  having  a volume  approximately 
5(X)0  times  larger.  The  experimental  results  reported 
herein  were  obtained  on  a heat  sink  which  was  3.75 
times  larger  Yolumetrically  than  the  original  mixlel. 
Tliis  report  has  presented  1)  the  experimental  results 
of  using  the  larger  model,  2)  a comparison  of  these  re- 
sults with  those  predicted  by  a previously  developed 
computer  program  and  3)  a comparison  of  results 
from  the  two  scale  model  test  series. 

Over  the  range  of  flow  rates  tested,  it  was  found 
that  a rather  well-defined  transition  developed  in  the 
pattern  of  water  flow  through  the  sink.  At  the  low- 
est flow  rates,  distinct  water  layers  formed  which  were 
essentially  isothermal  and  melting  progressed  primarily 
from  the  top  downward,  rather  than  radially  inward 
from  the  annular  space  around  the  ice  cylinder.  At 
the  higliest  flow  rates,  the  dominant  melting  mode  was 
radial  and  the  flow  is  rather  well-mixed.  At  the  inter- 
mediate flow  rates  the  melt  pattern  represents  a bal- 
ance between  mixed  and  stratified  flow. 

Except  for  some  anomalous  behavior  at  the  inter- 
mediate flow  rates,  the  average  outlet  water  tempera- 
ture during  the  melting  period  tended  to  increase  as 
flow  rate  increased.  For  both  the  model  studies  tliis 
average  value  vaned  between  about  36°F  and  44‘U-'; 
subsequent  to  melting  the  sink  outlet  temperature  in- 
creases rapidly.  The  coolant  water  flow  fate  did  not 
appear  to  significantly  influence  the  rate  of  melting, 
particularly  during  the  last  half  of  the  melting  period 
(Fig.  13).  During  the  early  stages  of  melting  the  influ- 
ence of  water  flow  rate  tended  to  be  masked  by  other 
factors,  such  as  a difference  in  the  average  ice  temp- 
eiature  at  the  start  of  the  test. 

A comparison  between  computer  predictions  and 
experimental  measurement  was  made  on  the  basis  of 

1)  mean  outlet  temperature  during  the  melting  period, 

2)  overall  time  required  to  complete  melting,  and  3) 
behavior  of  the  outlet  water  temperature  after  all  the 
ice  is  melted.  The  computer  solutions  for  outlet  water 
temperatures  during  average  melting  periods  tended  to 
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predict  a sliglitly  higlier  lenipcraturc  than  was  actu- 
ally measured  by  about  1°  to  in  most  cases  (Fig. 

1 1 ).  this  IS  I'elt  to  represent  a satisfactory  prediction, 
and  the  overestimation  of  temperature  tends  to  be 
sliglitly  conservative.  Also,  as  presently  structured, 
the  computer  program  fails  to  simulate  the  tempera- 
ture stratification  phenomenon  observed  subsequent 
to  melting  in  the  lowest  flow  rate  expenment.  Rather, 
the  program  predicts  a sliglit  plateau  in  outlet  temp- 
erature which  is  more  representative  of  the  higlier 
How  rale  tests.  Tliis  limitation  is  not  felt  to  be  a 
serious  drawback,  and  a simple  relationship  is  pre- 
sented in  this  report  to  more  accurately  depict  this 
temperature-time  relationship. 

It  appears  that  the  computer  predictions  could  be 
used  to  design  a prototype  sink  whose  actual  outlet 
water  temperature  will  be  generally  lower  than  that 
predicted.  Tlie  small  error  in  the  prediction  could 
be  considered  to  be  a sliglit  safety  factor  in  the  selec- 
tion and  si/.ing  of  a prototype  system.  The  computer 
program  did  predict  the  melting  limes  rather  accur- 
ately (Table  VII). 

Ttie  scaling  relationships  presented  in  Appendix  B 
appear  to  provide  a reasonable  methrxl  for  compar- 
ing the  results  ol  the  two  mixlcl  heat  sinks  of  differ- 
ent si/e.  At  present  there  appears  to  be  no  reason  to 
doubt  the  use  of  these  scaling  relationships  to  esti- 
mate possible  heat  rejection  and  coolant  water  flow 
rates  in  a prototype  installation.  In  those  flow  rate 
regions  where  good  correlation  between  the  models 
was  nut  obtained  using  these  scaling  relationships,  a 
iiKthod  was  developed  to  account  for  nonlinearities. 
h is  recommended  that  the  relationships  of  Appendix 
B should  be  employed  to  detemiine  initial  values  in 
tlie  sizing  of  a prototype  sink,  and  that  these  values 
should  be  refined  by  the  techniques  described  in  the 
section  of  this  report  which  discussed  the  approxima- 
tion of  the  heat  transfer  coefficient.  Based  on  the 
small  number  of  tests  conducted  during  these  pro- 
grams, it  is  believed  that  the  overall  results  indicate 
the  general  applicability  of  the  scaling  relationships 
to  approximate  the  behavior  of  any  size  heat  sink 
used  in  a similar  manner,  giving  due  consideration  to 
the  nonlinear  effects  indicated  by  the  Reynolds  num- 
ber analysis. 

By  using  a simplified  radial  melting  mathematical 
model,  variations  in  gross  average  heat  transfer  coef- 
ficients were  obtained.  When  the  heat  transfer  coef- 
ficient was  calculated  for  lO-hour  intervals,  the  results 


were  found  to  correlate  quite  well  lor  the  lower  flow 
rate  tests  with  values  determined  by  Vainer'*  who 
studied  natural  laminar  convective  melting  of  ice  ori- 
ented in  a vertical  plane.  Heat  transler  coefficients 
were  correlated  with  the  flow  conditions  by  using  a 
Nusselt  number  based  on  the  heiglit  ol  the  ice  and  a 
Reynolds  number  based  on  an  equivalent  diameter  of 
the  annular  space  between  the  ice  and  the  heat  sink 
tank.  Tlie  correlation  indicated  that  two  distinct  zones 
exist,  one  for  test  flow  rates  of  4.73  gpin  or  less  and 
the  other  for  flows  of  7.50  gprn  or  greater.  Tlie  transi- 
tion between  the  zones  appears  to  occur  when  the  Rey- 
nolds number  dunng  the  initial  hours  of  melting  lies  in 
the  range  of  KX)  to  120.  Tlie  two  relationships  devel- 
oped were 

Nu„  = 0.881(Re„f  Pr)*’  *'^^ 
for 

450  < (Reij^.Pr)<  1400 
and 

NUp  = 39.72(Re„^.  Pr)"  ^'** 
for 

1400  < (ReDKPr)<  3500. 

Tlie  Reynolds  number  analysis  was  also  found  to  be 
useful  in  explaining  the  nonlineanties  that  occurred 
when  using  the  flow  rate  scaling  prtKcdurcs  of  Appen- 
dix B. 

Within  reasonable  limits  of  accuracy,  the  computer 
program  does  simulate  the  behavior  of  the  heat  sink 
and  can  be  employed  as  a basis  of  design  for  a proto- 
type sink.  Tlie  use  of  the  computer  program  as  a de- 
sign basis,  rather  than  the  scaling  of  the  experimental 
results,  is  recommended  since  a prototype  installation 
would  employ  both  variable  coolant  flow  rates  and 
heal  rejection  rates  - conditions  that  are  easily  accom- 
modated using  the  computet  model. 

From  the  results  of  the  expenmenlal  tests  conduc- 
ted using  the  Model  I and  Model  II  heat  sinks,  it  is 
recommended  that  future  work  should  primarily  be 
concerned  with  the  investigation  of  alternative  inlet/ 
outlet  header  configurations  and  the  behavior  of  the 
sink  at  high  coolant  water  flow  rates.  Tlie  use  of  a 
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spray  type  inlet  header  in  a limited  number  of  Mrxlel 
I tests  resulted  in  the  creation  of  very  favorable  con- 
ditions and  could  result  in  an  improvement  in  the  ef- 
ficiency of  the  power  plant/heat  sink  system.  It  is 
felt  to  be  important  to  analyze  the  effects  of  this  and 
other  proposed  configurations  prior  to  the  construc- 
tion of  a prototype  sink. 

The  tendency  of  the  sink  to  provide  a constant  av- 
erage outlet  water  temperature  during  the  melting 
period  in  the  higlier  flow  rate  tests  also  appears  to 
merit  further  investigation.  Should  this  trend  contin- 
ue, extremely  simple  relationships  can  be  developed 
to  accurately  predict  the  behavior  of  the  sink  at  these 
higlier  flow  rate  levels. 

Wliile  mentioned  only  briefly  in  this  report,  the 
authors  would  like  to  emphasize  that  to  date  little 
information  can  be  found  in  the  literature  concerning 
tire  heat  and  mass  transfer  at  an  ice/water  interface 
under  conditions  similar  to  those  which  iKcur  in  the 
heat  sink. 
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Figure  A2.  Temperature- time  and  rate  of  melting  curves  for  the  2.50-gpm  test. 


100 


Time,  hrs 

a.  Inlet  and  outlet  water  temperatures. 


h.  Change  in  water  table. 


Tim«.  hr« 

c.  Rate  of  melting. 


Figure  A4.  Temperature-time  ami  rate  of  melting  curves  for  the  7.5(Fgpm  test. 
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Figure  A 5.  Temperature-time  and  rate  of  melting  curves  for  the  13. 50-gpm  test. 
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APPENDIX  B:  HEAT  SINK  SCALING  AND 
SIMILARITY  RELATIONSHIPS 


(Eroiii  Brown  and  Omnn’ ) 


Dimensional  scaling 

Tlic  dimensional  scaling  is  relatively  straiglitfor- 
watd.  Quite  simply  tor  the  prototype  and  model  to 
be  comparable : 


A )|nodel 

(Bl) 

~ if^/^a)proto 

(B2) 

where  = the  initial  radius  of  the  ice  cylinder 
R = the  radius  at  comparable  times 
/,  = the  initial  length  of  the  ice  cylinder. 

Temperature  scaling 

Meat  IS  transferred  to  the  ice  by  natural  convection 
which,  particularly  for  water,  is  higlily  dependent  on 
both  temperature  difference  and  absolute  tempera- 
lure.  Tins  dependence  IS  due  to  the  density-tempera- 
ture relationship  of  water  which  exhibits  a maximum 
value  at  4°C.  Since  the  heal  transfer  pr(x.ess  is  so 
higlily  dependent  on  temperature,  the  temperatures 
of  the  rmxlel  and  the  prototype  must  be  directly 
equivalent.  Tlius. 

T =T 

'model  proto 


Tins  allows  for  the  formulation  of  a “heat  storage" 
parameter; 


/V  =- 


Qh 


And  for  the  prototype  and  model  to  be  comparable 


^^s^model  ” ' 


We  have  a requirement  that  the  heat  rejection  rates 
must  be  related.  Tints,  the  percentage  of  heat  rejected 
tlirougli  melting  must  be  comparable  for  the  model 
and  prototype.  Tlius, 


or 


^ m(nn>del)  _ ^R(inodcl) 
^m(prolo)  ^R(proU>) 


and 

•^^mirdel  “ ■^^prulo- 

Heat  transfer  scaling 

For  the  heal  transfer  phenomena  in  the  mixlel  to 
be  comparable  to  those  of  the  prototype,  the  relative 
effects  of  the  vanous  modes  of  heat  transfer  and  stor- 
age must  be  equaled.  In  this  particular  study,  heat  is 
absorbed  in  the  heal  sink  by  two  methods:  1 ) heal- 
ing of  the  water  m the  sink  , and  2 ) melting  ol 
the  ice.  thereby  utilizing  its  latent  heat 

Tlius,  the  total  heat  rejected  to  the  sink  up  to 
a given  point  in  lime  is 


Now 

= ''conv^<7'B-’2) 

where  = convective  heat  transfer  cocificient 
A = surface  area  of  ice 
Tj,  = bulk  water  temperature. 

Tlius 

^^erwrv^  ^2*)nnidel 

^R(pn>l(>) 

As  noted  above,  the  temperatures  must  remain  equiva- 
lent, thus 
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^R(miHld)  _ ^^*atnv  ^model  _ ^^*cf»nv  ^0  ^ntKJd 
^ R(pn)t<>)  ^^^conv  V^olo  ^^^lonv  ^0  V^oto 

Fi>r  turbulent  !l(»w.  is  iiulcpoiidciil  til  si/c;  thus 


^mculel  Vf'dn 


^^citnv(m(Kk’l)  ^^c(»nv(prt»l(»)' 


Tliese  scaling  rclalitinships  can  now  be  sumniari/cd. 
For  the  general  case: 

^R(m(Kiel)  _ ^^^c«nv  ^modcl  (B3) 

^R(pnilc>)  ^^*cimv 

For  the  turbulent  case: 

^R(m<>del)  _ (B4) 

^R(pntl<>)  "^pnito 

where  A - ItiRq!..  A means  lor  scaling  the  time  fac- 
tor must  ntiw  be  derived.  For  this  case,  tlic  percentage 
of  total  heat  abstubed  must  be  equal  for  equivalently 
scaled  times.  Tlius 


Cr  ^ r Qjo 
Qjc  J eTc 


wliere  C^tc  is  the  total  heat-absorbing  capacity  of  the 
heat  sink  system. 

Assuming  that  heat  is  rejected  at  a constant  rate, 
die  coolant  water  flows  at  constant  volume,  and  the 
initial  sink  consists  solely  of  ice,  the  following  rela- 
tionship holds: 


Cb  0 AT 

— =N=  2 

Qj,.  ^ (//  + C,Ar„„) 


where  AT  = constant  temperature  difference  across 
the  sink,  r„-r„„, 

= maximum  allowable  sink  temperature, 
T -32°F 

W = mass  flow  rate,  Ibm/hr 
Cp  = specific  heal  of  water,  1.0  Btu/lbm  °F 
Af^i,  = mass,  lb 

//=  latent  heat  of  fusion,  144  Btu/lbm 
0 = lime,  hr. 


we  may  obtain 


■(— ) 
v^sink  'prolo 


\ ''^sink /model  V 'sink/prolo 

We  note  that  this  is  the  same  result  as  would  be  ob- 
tained by  considering  that  the  percentage  of  total 
mass  circulated  througli  the  sink  at  equivalent  times 
must  be  comparable  Ttius 

r 0 1 =r_o 1 (B6) 

Kink/H  model  Kink/H  proto 

It  should  be  noted  that  this  result  is  also  obtainable 
from  the  “storage"  parameter: 


where  = (2itrL)pHr 

p = density  of  ice,  Ibm/fl^ 
r = drIdO. 


{2nrL)pH‘^ 
c/^ 

WCy,AT 


For  the  model  and  prototype  to  be  considered  equiva- 
lent 


N = N 

' s(model)  ^*s(prolo)  ‘ 


Tlius,  eliminating  equivalent  terms. 


'I.RR\  ^ /l.Ri 

. W /model  \ W 


where  R is  the  instantaneous  radius  at  time  d.  Now 
consider 


For  comparable  cases  (of  model  vs  prototype)  the 
temperature  AT  and  are  equal.  Thus,  since 


It  is  essential  that  a dimensionless  time  factor  be 
formed  if  equivalent  time  scaling  is  to  exist: 
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\0*/mcMlt’l  w*/pr<»lo 


where  W*  is  lime  in  wtiicli  all  ice  will  be  melted  and 
sink  water  is  broiiglit  to  its  maximnm  allowable  temp- 
erature . and 

Oh 


<7(0/0*) 

where  /(y  is  the  initial  ice  cylinder  radius.  Now 
<«''«o)n,..dcl  = («/«o)pro.o 

when 

<^'/«*)n,..det  = f«/0*)pr,.,<. 

Til  US 

^ r<7(/7//(o)] 

|_</((;/0*Jm..dcl  l_(/(0/0*  'Jprnto 

This  result  is  used  in  obtaining 

I r / /( 

k'  '\R„  '7\  0*1  </(0/0*)J model 

k'[\«o  7\0*/  <7(o7o*ljprolo 
Canceling  equivalent  temis 

\kW*/mi>del  vt^* 'prolo 


Summary  of  dimensional  and  heal  transfer 
similarities 

Based  upon  the  above  discussion  when  the  follow- 
ing series  of  requirements  are  followed,  the  scale 
model  and  prototype  are  comparable: 

I.  Temperature  scaling. 


•^^model  •^^prolo 


T = r 

model  prolo 


2.  Dimensional  scaling: 


W«0)model  = (/^/«0  Vfoto 


(Rq/I-  Imodel  ” ^Rplf-  ^proio 


3.  Time  scaling: 


(0/0*)modd=<0/»*»pro.o 

4.  Ileal  rejection  rate  scaling: 

^ R(niodel)  ^^^conv  ^ ^odel 
^K(prolo)  ^^^conv  ^ ^roto 


where  A is  Zirr^L. 


Development  of  the  scaling  relationships 

extrapolation  to  the  prototype  may  be  accom- 
plished by  plotting  7j„  and  vs  (010*)  and  (r/r^) 
vs  (010*):  these  relationships  are  applicable  for  a 
specific  ATfri^-  7]^,)  condition.  Tlie  lime  scale  is 
established  by  selection  of  the  flow  rale  (WjM): 


which  becomt's  the  same  result  as  obtained  above 
wlien  inverted  and  multiplied  by  trp^.^: 


y M /miMlel  \ M /pi 


® (model)  _ \^R  /model 
®*(prol<.)  MinkX 

\^R  /prolo 


("+f’p7m«) 


(//+r  T ) 

' p ‘m»x  ' 
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Ni)w 

+ ^ r '/nK«  >n..,del  = ^max  >pr..t<. 


S(^ 


®*(mmk-l)  _ (Pi"^o^  /^R  )m<>del 
^ (prcdo)  {p(nRlLl. 

R ^proto 


_ (^O^)model  ^R(prol<>) 

(^0^'Vfolo  ^R(niodel) 

since 

^R(prnt(>)  _ ^^*c(>nv \>r(>lo  _ ^^conv  ^ ^^0^' V*‘t>**> 
^^R(imidel)  ^^’conV^^modcl  ^^conv  ^^^O^'^model 


Tlie  following  relationship  holds  between  model  and 
prototype  for  the  initial  sink  radius,  relative  to  tlie 
total  useful  sink  life: 


0* 

(model) 

0* 

" (prulii) 


_ (''o^-^model  . 

^^0^'  VnUo 

_ ^O(miidel)  ^ 
^Otprolo) 


^^*conv  ^0  )pruli) 
(^'conv^'o^-^model 


conv(prol(>) 


conv(modei) 


(B8) 


A.S  noted  previously  for  turbulent  How, is  inde- 
pendent of  siire;  thus 


^*conv(pr(>lt>)  “ ^*conv(m<»deI)' 

It  should  be  recognized  that  for  a scale-model  test 
expenencing  laminar  flow,  any  extrapolation  to  a 
large  size  sink  undergoing  turbulent  heat  transfer  will 
residt  in  conservative  estimates  of  coolant  water  heat- 
ing rates.  Turbulent  How  will  induce  greater  mixing 
and  will  thus  produce  higtier  heat  transfer  coeffici- 
ents (htglier  Nusselt  numbers).  Consequently,  ex- 
penmentaJ  results  obtained  under  laminar  flow  will. 


wtien  extrapolated  to  a large  sink,  tend  to  err  on  the 
conservative  side  under  the  similarity  relations  de- 
rived above. 

Tlius 


^ (model)  _ ^0(m(»del) 

n*  Q 

(proto)  ^0(pr(»to) 


and 

h = ^^0^' Vr»‘»  n 

^R(proto)  ip  w . ^R(model) 

(Aq''  ^model 


III  SO 


WO 


A/. 


sink/model 


WO* 


sink/proto 


W. 


proto 


A/.; 


sink/mode! 


M. 


sink 


w^ 


model 


/proto 


Rr 


sink  'model 
^0 


M. 


'sink 


W, 


model 


0 'proto 


nRll. 


nRl/. 


nn»del 


0‘'/nit>del 


0 proto 


proto 


(/?0 1'  )ni<)del 


'^model 


(Bd) 


Using  the  above  relationships,  the  dimensional  and 
heal  transfer  similarities  are  satisfied.  For  the  series 
of  experiments  conducted  during  the  study  of  the  ice 
water  heat  sink  concept.  Table  Al  provides  a compari- 
son between  the  scale  models  and  the  prototype  based 
on  nominal  tank  dimensions. 


Table  Al.  Comparison  between  the  scale  models  and  the  prototype 
based  on  nominal  tank  dimensions 


Pnttotype  sink 
fftSftiJiam  X HO  ft) 

MtnJel  11  sink 
(f>  ft  iham  X to  ft) 

Model  1 sink 
(4  ft  diam  X 6 

Heat  rciectiori 

5.96X  10* 

47,763* 

19.105 

rale  (Htu/hr) 

5.22X  10* 

42.013 

I6,H0S 

2.61  X 10* 

21,006* 

H,405 

Cooliint  water 

1,6K0 

13.50 

5.40 

How  rate  (gpm) 

I.24S 

10.00 

4.00 

935 

7.50 

3.00 

590 

4.73 

1.R9 

310 

2.50 

1.00 

205 

1.55* 

0.62 

* No  tests  vb«re  conducted  at  this  value. 
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\PPFNI)IX  C DERIVATION  OF  THE  RELATION 
SHIPS  FOR  THE  HEAT  TRANSFER  COEFFICIENT. 
REYNOLDS  AND  NLSSELT  NUMBERS 


Rckiiiimships  are  developed  here  to  approximate 
the  heat  transier  eoelTieient.  Reynolds  number  and 
Nusselt  number.  Sample  calculations  are  also  included 
to  illustrate  the  application  of  these  relationships. 

Heal  transfer  coefficient 

A first  order  approximation  of  the  coefficient  of 
heal  liansfer  can  be  made  by  considering  the  gross 
melting  priKCss  which  occurs  within  the  heal  sink. 
Considering  only  that  portion  of  the  total  heal  re- 
jected to  the  ice  dunng  any  finite  period  of 
lime,  we  have 

JQ^  = hA  ST  di 

where  li  = heat  transfer  coefficient 
A = area  of  the  melting  surface 
= temperature  difference  between  the  melt- 
ing surface  and  the  bulk  water  tcmjierature 
Jl  = length  of  the  time  period. 

If  the  ice  is  initially  at  32°F  then  all  the  heal  re- 
jected must  be  used  to  melt  the  ice.  If  it  is  assumed 
that  only  radial  melting  (xrciirs  and  melting  is  uniform 
along  the  height  of  the  ice  cylinder,  then 

= ! ‘IM, 

* m I 

where  Af,  = p,rr/?^//  Tliercfore 

= pj.  IvRHdR 

where  Pj  = density  of  the  ice 

/-  = latent  heat  of  fusion 

R = radius  of  tlie  ice  cylinder 

//  = lieiglit  of  the  ice  cylinder  (constant). 

-Since  melting  is  assumed  to  occur  only  on  the  cir- 
cumferential surface  of  the  ice,  the  factor  of  InRH  is 
tlie  effective  melting  surface  .irea.  Therefore  the 
equation  becomes 


dQ„^=pJ.AdR. 

From  this  expression  and  the  relation  involving  the 
heat  transfer  coefficient,  the  following  expression  is 
obtained. 

(cri 

ir„ 

Tims,  the  heat  transfer  ctxrfficient  can  be  approxi- 
mated when  the  rale  of  radial  melting  and  the  bulk 
water  temperature  are  known.  Since  both  these  para- 
meters vary  with  time  and  elevation  in  tlie  sink,  a spe- 
cific detemiination  of  the  heat  transfer  coefficient  is 
not  possible. 

To  approximate  the  heat  transfer  coefficient  it  can 
be  assumed,  as  a first  order  approximation,  that  the 
bulk  water  temperature  would  be  equal  to  the  average 
value  of  the  inlet  and  outlet  water  during  the  melting 
|x:riod.  Tliat  is 

(C2) 

where  Tg  = average  bulk  water  temperature  during 
melting 

= average  inlet  water  temperature  during 
melting 

= average  outlet  water  temperature  during 
melting. 

From  this  the  mean  temperature  difference  between 
tlte  ice  and  the  water  is 

ST^  = t^-n.  (C3) 

Considering  the  melting  prrKess  with  respect  to  a 
time  interval  of  duration  St  hours  the  rate  of  change 
in  radius  can  be  approximated  by 

dR  ^ SR 
dt  St 
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wliorc  AK  IS  ilio  (.liaiige  in  ihc  liulnis  (.liirini;  the  tune 
interval  A/.  Tlnis  tlic  approximation  lot  the  heafeo- 
ettieient  beci'ines 


h 


A/v’ 

('(r  +/■ 

I ^ ' in  out ' I 


((•4) 


toi) 

4 ‘ ' 

where  lt'=  mass  Ilow  rate  (Ihm/hr)  aiul  - ilensity  ol 
water  (IHin/ft’ ). 


llie  appheatioti  ol  tins  ecpiation  is  not  restiieted  to 
any  s|K'eil'ie  time  inters al  ihiiation.  It  Ar  is  selceteil  to 
he  the  total  meltint;  period,  then  \R  would  he  the 
ehange  m ladms  Iroiii  the  mitud  value  to  zero.  It  Ai 
is  less  thati  the  total  melting  penod,  the  change  in  the 
radius  lor  that  penod  must  he  approximated  using  ad- 
ditional information,  such  as  the  rate  of  melting. 

Reynolds  numher 

Nondimcnsional  parameters  have  heen  employed  to 
establish  a generalized  correlation  between  the  heat 
transfer  ciK'lTieient  and  the  coolant  water  Ilow  rates. 
Water  velocity  effects  arc  normally  considered  by  use 
of  a Reynolds  number  Re  where 


where  U = mean  water  velocity  (ft/hr) 

D=  characteristic  dimcnsioti  (diameter)  (It ) 
u = kiticmatic  viscosity  (It^/lir). 

For  flow  in  annular  passages  the  characteristic  dim- 
ension norttially  selected  is  the  “cc|uivalent”  diameter 
of  the  annul'is"'  which  is  defined  by 


4^(Dl-Dh 

Dl  = 

rt(/^2  *D, ) 


(C.S) 


Subsliluting  these  eipiations  into  the  expiession  lor 
the  Reynolds  number  yields 


4 If 

t^P^nt/Ji-t-/;,) 


((  7| 


where  /),  is  nominal  lank  diamelei  (It ).  and  /)^  is  the 
average  ice  diamelei  during  the  lime  inleixal  ol  a|'pio\i- 
malion. 

Tlie  Reynolds  numher  is  appio.ximaled  loi  the  same 
series  of  time  intervals  as  l!ie  heal  Iranslei  coellicieni 


Nussell  number 

Tlie  Nussell  niimbei  is  employed  to  noiulimension- 
alize  the  heal  liaiisfer  coellicieni  loi  coiiel.iiion  with 
the  Reynolds  luimbei . This  lelalionship  is 

/r//„ 

Nu,  = -^  (CK) 

where  /i  = heal  transfer  coeflicienl  (Ulii/hi  ir  “l  ) 

//y  = initial  heiglit  of  the  ice  cy  hndci  (It ) 

A = ihermal  conductivity  ol  the  water  (Ulii/hr 
ft  °F). 

Tlie  use  of  the  initial  ice  cylindei  height  as  the  chai- 
acleristic  dimension  agrees  with  the  assumptions  em- 
ployed in  the  approximation  of  the  heat  liansler  cik'I- 
licienl  and  appears  to  provide  a leasonable  coiielalion 
Tlie  basis  for  this  selection  was  discussed  pieviousls 


where  />/:  = equivalent  diameter  ( ft ) 

= Ollier  diameter  o(  the  annulus  (It) 

/.>!  = inner  diameter  ol  the  annulus  (ft). 

In  the  ice-water  heat  sink,  Dj  is  the  diameter  of  the 
lest  apparatus  (constant)  while  O,  is  the  average  di- 
ameter of  the  ice  cylinder  during  the  Interval  of  ap- 
proximation (vanable).  Tlie  mean  water  vehxiity  Is 
determined  by  dividing  the  volumetric  Ilow  rale  by 
the  average  cross-sectional  area  of  the  annulus.  Tlius, 
the  mean  flow  velocity  is  calculated  from 


fsample  calculation 

To  illustrate  the  application  ol  these  lelalionships, 
the  heal  transler  cociricienl,  Reynolds  and  Nussell 
numbers  for  the  first  lO-hour  interval  of  test  IH  w ill 
be  calculated.  Tlie  data  tised  in  these  sample  calcula- 
tions can  be  foiiiid  in  Table  II. 

//car  transjcrcocjjicicnt.  From  equation  ("4 


/I  = 


P|/. 


^ 

1(7-  +7'  ,)-12 

' 111  ' nut  ' 


A/ 
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Using  die  constants 


Pj  = 56.7  Ibm/ft’ 

/.  = 144  Btu/lbin 

and  the  following  data  for  the  firs'  10  hours  of  test  1 B 
7-,n  = 54.2°F 


algebraic  average  of  the  bulk  water  temperature  and 
the  temperature  of  the  heat  transfer  surface  (i.e.,  the 
ice).  For  diis  example  the  mean  fluid  temperature  7} 
is 


7}  = 


J_ 

2 


(46.0+32.0)°  F 


39°  F. 


37.7°F 


AR  = 0.34  ft 


Ar=  10  hr 


Tltc  thermal  conductivity  is  evaluated  at  39°F,  so  that 
die  Nusselt  number  is 

_ (19.90  Btu/hr  ft^  °FK8.68  ft) 

““  (0.331  Btu/hr  ft^  °F) 


(56.7X144)  0.34 

i(54.2+37.7)-32 


Btu/hr  ft'  °F 


or 


Nup  = 522. 


or 


h=  19.90  Btu/hr  ft'  °F. 

Reynolds  number.  Tlie  nominal  test  (low  rate  was 
4.73  gpm  or  37.94  ft'/hr.  During  the  first  10  hours 
the  ice  radius  changes  from  3.00  ft  to  2.66  ft  (AR  = 
0.34  ft)  so  that  die  average  ice  radius  is  2.83  ft.  Tlie 
kinematic  viscosity  is  evaluated  at  the  bulk  water 
temperature  of  46.0°F,  thus,  the  following  constants 
apply: 


37.94p^  (Ibm/hr) 

D,  = 6.00  (ft) 

Di  = i(6.00+5.32)=  5.66  (ft) 

u = 0.0545  (ft'/hr). 

Thus  from  eq  C7 

_ (4X37.94  p^) 

' irp^ (0.0545X6.00+ 5.66) 

Rcpj  « 76.0 


Nusseh  number.  The  height  employed  in  eq  C8  is 
8.68  ft  and  the  thermal  conductivity  of  the  water  is 
evaluated  at  the  mean  fluid  temperature,  which  is  the 
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M’l’I  NDIX  1)  DFTtRMINATION  OF  FREEZING 
R \ri  S \M)  REFRIGERATION  LOADS* 

Hii'  luwiiii;  uci'iiii’lty  of  llio  heal  sink  is  illustratcil 
ni  I i.'iiK  l»l  Flic  use  ol  the  ciHilinp  coils  along  the 
lank  ssall  i.snlis  in  railial  Iree/ing  inward  front  the 
wall.  I ree/ing  at  the  holloin  ol  the  lank  and  end  cl- 
tecis  at  the  U|'per  sniface  have  been  neglected  in 
Older  to  develop  a siinphUerl  prevliclion  lechnic|ne. 

Pie  water  in  the  tank  is  assumed  to  be  initially  at 

With  these  assumptions  the  I'lec/'ing  of  the  sink  is 
dependent  only  upon  the  ladial  coordinate,  and  thus 
the  lempeialnre  distribution  m the  ice  annulus  is 
governed  by  the  following  partial  differential  ei|iia- 
lion 


llie  associated  bonndarv  conditions  arc 


/(/f,,/)=  To 


nR. 


where  R^^  = radius  at  ice/water  interface.  Tlie  heal 
transfer  at  the  ice/waler  interface  is  described  by  the 
following  relationship 

3/  p,/.  dR  K 

Since  the  variation  in  ice  temperature  with  time  is 
slow,  a steady  state  solution  to  the  differentitt]  equa- 
IKin  may  he  used.  Tlius  the  temperature  at  any  time 
t IS  given  by 

(7-|-7-„)ln(^) 

•r(r)=7„+— - 

In/ 


\"  7 


(Dl) 


Upon  substitution  of  this  relationship  into  the  heat 

transfer  relationsfiip,  the  following  expression  for  the 

time  reiiuired  to  form  an  ice  annulus  of  radius  R 
' m 

results: 


-2  Ini 


(D2) 


• Ortfunal  Jenvation  of  cquatians  hy  R.  (irmde,  ref- 
erenev  7. 


For  liiC  Model  II  heat  sink  the  following  aic  the 
parameter  values  to  be  employed 
= .VO  ft 

Pl  = .S6.7  Ihm/ft^ 

/.  = 144  Btu/lhm 


7-  =-70°F 
^ = 1.26  Btu/hr  ft  °f'. 

Tliese  values  were  substituted  into  the  above  equa- 
tion and  Figure  D2  was  developed  to  show  the  rate  of 
growth  of  the  ice  annulus  with  time.  Tliis  relationship 
predicts  that  free/ingof  the  Model  II  sink  will  iKCur 
within  6 days,  which  agrees  with  the  measured  free?e- 
up  time.  Althougli  the  tank  wall  temperature  was  not 
-70°F  at  the  start  id' free/.ing,  only  3 to  4 hours  was 
required  to  attain  that  temperature. 

It  should  also  be  noted  that  the  water  witltin  tlie 
tank  was  agitated  hy  the  air  bubbling  system  and  that, 
altliougli  the  equation  used  to  estimate  freezing  times 
IS  for  freezing  along  a quiescent  boundary  , it  predicted 
the  overall  freezing  lime  rather  well. 

An  estimation  of  the  rate  ol  heat  flow  to  the  refrig- 
eration coiis  can  be  made  from 


Q=2vRJI„k 


<iR 


H 


which  upon  substitution  ol  ihe  expression  tor  T 
becomes 


c>  = 


2tr//oA(7-,  - To) 


LIsing  a value  of = 9.0  ft  (a  typical  test  dimen- 
sion) Ihe  variation  in  heal  (low  with  respect  to  time 
has  also  been  plotted  in  Figure  D2.  After  the  First 
day  the  rate  of  heat  flow  has  decreased  from  more 
tlian  68,000  Btu/hr  to  approximately  20,000  Btu/ 
hr  as  a result  of  formation  of  a 1-ft-thick  annular 
section  of  ice. 


I 
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JWCKOIWO  PAOI  fiUMUNOT  nu«D 


It  hatl  been  »nj;ii\all>  assuincil  that  tbe  length  of 
the  freezing  [>eriml  eoiikl  be  substantially  reduced  by 
preloading  the  sink  with  ice  blocks. however,  a nioie 
detailed  examiiiaiion  tended  toconllict  with  this  Itne 
of  reasoning  An  exairiinalioii  of  l igurc  1)3,  wliich 
sliows  the  theoretical  rale  of  freezing  based  on  vol- 
ume, indicated  that  .SI9f  of  the  mass  of  the  Model  II 
sink  wa.s  already  liozen  at  the  end  of  the  first  day. 

It  should  be  noted  that  when  the  sink  was  preloaded 
with  approximately  50'^  ice.  the  length  ol  the  trecz.iiig 
period  was  reduced  by  only  one  day.  Tlius.  the  ob- 
served decrease  coiielatcs  well  with  the  predicted  be- 
havior, even  ihougti  two  “distinct"  sets  of  freezing 
conditions  exist.  I iiriher  examination  of  Figure  1)3 
indicates  why  both  of  these  zones  agree.  Note  that 
during  the  last  day  of  (recz.ing  only  3%  ol  the  total 
volume  becomes  frozen;  during  the  last  two  days 
imly  \07r  of  the  total  volume  can  be  frozen,  even 
thougfi  a wall  Icinperature  of -70°  is  used.  This  is  a 
result  of  the  long  path  througli  which  the  heal  must 
he  conducted,  which  during  the  last  two  days  varies 


liom  2.1  to  3.0  ft  at  the  completion  ol  freezing.  Iri- 
dcpendcnl  of  whether  or  not  the  outer  layer  was 
composed  of  ice  hlwks  or  frozen  water,  the  length 
of  this  path  rematiis  the  same.  Tlius  a substantial 
decrease  in  the  length  of  the  lieezmg  period  is  pos- 
sible only  if  the  length  ol  the  heal  conduction  path 
Ls  made  shorter  or  a massive  initial  ice  load  is  used. 

To  decrease  the  ireezeup  period  by  507f.  it  would  be 
necessary  to  have  an  iiiitiaJ  load  of  ice  equal  to  827? 
ol  the  total  mass  of  the  sink.  It  should  be  noted  that 
in  previous  solid  cylinder  tests  the  peri.cntagcof  ice 
mavs  was  on  the  order  ol  S.ST?  to  '■HYX  ol  total  mass, 
and  thus  pieloauing  blocks  to  achieve  an  827?  load 
factor  would  be  impossible. 

Tins  IS  not  to  say  that  preloading  with  ice  is  with- 
out benefit.  In  Figure  1)2  it  can  be  seen  that  the 
greatest  heat  load  on  the  refrigeration  system  occurs 
dunng  the  same  period  that  the  miiiaJ  50'??  of  the 
sink  mass  is  Irozen.  Tims  preloadmg  of  ice  can  be 
translated  into  reduced  initial  heat  loads  and  reduced 
refrigeration  costs. 
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APPENDIX  E:  APPROXIMATION  OF  STRESSES 
IN  THE  HEAT  SINK  TANK 


In  the  previims  ice-watcr  heat  sink  study it  had 
been  observed  that  fractures  occurred  in  the  ice  after 
freezing  had  (teen  completed  and  the  heating  cycle 
started.  IXiring  assembly  of  the  Model  II  lank,  bond- 
ed strain  gages  were  installed  to  monitor  induced  wall 
stresses  during  the  thawing  mode.  However,  because 
of  damage  during  tank  assembly  only  one  measuring 
gage  and  the  temperature  compensating  gage  were 
operational  during  the  lest  program.  Tlte  temperature 
compensating  gage  was  mounted  on  a section  of  steel 
of  the  same  type  as  used  to  make  the  heat  sink  tank 
and  exposed  to  the  same  temperature  conditions  as 
the  tank  wall.  When  connected  in  a bridge  network 
with  the  measuring  gage,  the  temperature  compensat- 
ing gage  eliminates  the  apparent  strain  readings  that 
result  from  temperature-induced  expansion  or  con- 
traction of  the  lank  wall.  The  active  gage  was  located 
at  the  midheight  of  the  tank  and  was  mounted  in  the 
circumferential  direction  to  measure  induced  hoop 
stresses.  A thermrxouple  mounted  at  the  same  loca- 
bon  on  the  inside  of  the  lank  was  used  to  measure 
tank  wall  temperature. 

As  an  approximation  of  the  true  stress  condition 
it  can  be  assumed  that  the  heat  sink  tank  experiences 
a stress  condition  similar  to  that  of  a long,  thin-walled 
pressure  vessel  without  edge  effects,  subjected  to  a 
uniform  internal  pressure.  Actually  the  Model  II  tank 
is  a rather  short  vessel  with  a length  to  diameter  ratio 
of  only  1.33  to  1,  but  since  the  gages  were  mounted 
at  midheight,  this  assumed  stress  condition  represents 
a reasonable  first-order  approximation. 

Under  such  a stress  condition  the  stress  in  the  wall 
can  be  expressed  as 


where  Og  = circumferential  wall  stress  (hoop  stress) 
P-  Internal  pressure 
R = radius  of  the  cylinder  (36.0  in.) 
t = wall  thickness  (0.375  in.). 


Since  from  Hooke’s  law 


o„  = hc„ 


where  /-.'is  Young’s  modulus  (3()x  10*  Ibf/in.^  for 
steel)  and  Cg  is  hoop  strain,  then 


(E2) 


To  establish  a base  value  for  strain  measurements,  a 
reading  of  the  strain  gages  was  taken  with  a water  head 
of  9 ft  and  at  a wall  temperature  of  44°F  (simulated 
operation).  Using  the  above  relationship  and  the 
strain  readings  during  the  heating  period  just  after  Ihe 
completion  of  freezing,  the  internal  pressures  as  a 
function  of  tank  wall  temperature  were  calculated. 

The  results  are  presented  in  Figure  F I . Also  shown  is 
the  variation  in  hoop  stress  in  the  tank  wall. 

Tlie  figure  indicates  that,  as  the  tank  wall  tempera- 
ture increases  from  -70°F  to  20°F.  there  are  signifi- 
cant variations  in  Ihe  internal  pressure  and  thus  the 
stress  induced  m the  wall.  As  its  temperature  increases, 
the  ice  cylinder  attempts  to  expand  in  volume;  how- 
ever, since  Ihe  tank  wall  limits  the  total  amount  of 
volumetric  expansion  possible,  significant  stresses  are 
induced  in  both  Ihe  ice  and  lank  wall.  Tlie  sudden  de- 
creases in  stress  which  are  evident  on  Figure  El  most 
likely  result  from  the  propagation  of  cracks  within  the 
ice.  Of  course  ice  cracking  is  not  restricted  to  the 
period  following  the  initial  heating  of  Ihe  ice;  localized 
stress  concentrations  can  produce  cracks  at  any  time 
during  the  heating  process. 

Wliat  is  worth  noting  in  this  experiment  is  the  mag- 
nitude of  the  induced  stresses.  Tlie  measurements  in- 
dicate maximum,  minimum,  and  average  internal  pres- 
sures of  100.4,  .34.3,  and  79.7  Ib/in.’,  respectively, 
and  corresponding  maximum,  minimum  and  average 
hoop  stres-scs  of  10,200,  3,300  and  7,620  Ib/in.*, 
respectively.  Of  course  these  stresses  and  pressures 
are  functions  of  the  experimental  lank  conditions, 
i.e.,  a 9-ft  water  head  and  an  initial  wall  temperature 
of44°F. 
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may  be  induced  during  the  lime  that  ilie  ice  is  hemp 
initially  healed  to  form  the  annular  gap  bciwcen  ibt 
tank  and  ice  cylinder  walls. 


Tlicsc  obseived  wall  stresses  are  significantly 
greater  than  ilie  nominal  unconfined  fracture 
strength  of  ice.  In  general  this  indicates  that,  in  the 
operation  of  a prototype,  significant  lank  stresses 


Figure  FI.  Variations  in  internal  pressure  and  tank  wall  stress  during  the  ice  heating 
mode. 
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